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Abstract 



Chiral lagrangians as effective field theories of QCD are successfully ap- 
plied to meson physics at low energies in the framework of chiral perturbation 
theory Because of their nonlinear structure these lagrangians allow for static 
soliton solutions which may be interpreted as baryons. Their semiclassical 
quantization, which provides the leading order in an 1/Nq expansion (Nc is 
the number of colors) turned out to be insufficient in many cases to obtain 
good agreement with empirical baryon observables. However with Nc = 3, 
large corrections are expected in the next-to-leading order which is carried by 
pionic fluctuations around the soliton background. The calculation of these 
corrections requires renormalization to 1-loop of the underlying field theory. 
We present a procedure to calculate the 1-loop contributions for a variety of 
baryonic observables. In contrast to chiral perturbation theory, terms with 
an arbitrary number of gradients may in principle contribute and the re- 
striction to low chiral orders can only be justified by the investigation of the 
scale independence of the results. The results generally give the right sign 
and magnitude to reduce the discrepancy between theory and experiment 
with one exception: the axial quantities. These suffer from the fact that the 
underlying current algebra mixes different Nc orders, which suggests a large 
and positive next-to-next-to-leading order contribution, which is probably 
sufficient to close the gap to experiment. 
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Chapter 1 
Introduction 



1.1 Preliminaries 

During the past 15 years, effective field theories have been proven to be 
valuable tools in different branches of physics. In hadron physics, the need 
for effective field theories is brought about by the intractability of QCD at 
low energies around lGeV. Of course, the distinction between QCD on the 
one hand and effective field theories on the other is somewhat artificial at 
the conceptual level since QCD itself is probably an effective field theory 
of something more fundamental. However, for practical purposes, it is clear 
that the examples of effective field theories we are considering, namely chi- 
ral lagrangians formulated in terms of pseudoscalar fields, are much more 
restricted in their domain of validity and richness of phenomena they can 
describe than QCD. 

The most prominent, i. e. energetically lowest lying, degrees of freedom of 
such hadron effective field theories are pions (and kaons) thought to be the 
goldstone modes of the spontaneously broken chiral symmetry of QCD. What 
makes these mesonic effective field theories tractable in a concise fashion is 
the concept of chiral perturbation theory (ChPT), pioneered by Weinberg 
[[I -10 and for the first time exploited by Gasser and Leutwyler || ||, which 
has been successfully applied to a variety of meson- and meson-baryon phe- 
nomena (for the latter case, see e.g. 0). ChPT allows to calculate processes 
in a way that starting from a meson lagrangian L systematically includes all 
relevant terms up to some power of external momenta, which correspond to 
the number of derivatives (chiral order) contained in the vertices generated 
from C, and the quark masses. The only theoretical input are universally 
accepted principles like chiral symmetry (and a certain pattern of symme- 
try breaking), Lorentz covariance etc. . The essential attraction of ChPT 
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stems from the fact that calculations are no longer model dependent in the 
traditional sense. Unless some of the underlying general principles could be 
proven wrong, and as long as higher degrees of freedom (e.g. vector mesons) 
can safely assumed to be frozen, one only has to decide on the power in ex- 
ternal momenta to which the calculation will be carried, and in this way the 
error introduced is under control. 

The price to pay for this conceptual stringency rests with the fact that even 
at small orders, there is a vast number of independent terms in the effective 
lagrangian, each comes along with its own low energy constant (LEC) to be 
fixed from experimental information. Since experimental knowledge, as well 
as the practical ability to include terms in a calculation are both limited, so 
is the applicability of ChPT. 

Baryons have been introduced into the overall framework of effective field 
theories in several different ways: One possibility, directly related to ChPT, is 
to treat them as external fields coupled in the most general chiral invariant 
manner to pions. A similar momentum (and quark mass) power counting 
like the one outlined above may then be employed to calculate the 'dressing' 
of the originally pointlike particle by meson loops, provided the baryon is 
treated as a heavy field (HBChPT) H. Unfortunately, in this form, the 
proliferation of terms allowed from chiral symmetry alone already at the 
lowest nontrivial chiral orders is even more worrisome than in the purely 
mesonic case, and consequently the applicability of HBChPT is restricted 
even further compared to meson ChPT. is lost. 

Therefore, the idea to create the baryons as solitons from the meson fields 
themselves without explicit fermionic degrees of freedom has been consid- 
ered an especially attractive alternative, which has the advantage that one 
might directly build upon the effective lagrangian established in ChPT in the 
mesonic sector. It is this aspect that defines the scope of our paper, namely, 
to investigate the extent to which the concept of baryons as solitons in ef- 
fective meson theories can be linked to the methods and results of purely 
mesonic ChPT. Evidently, the consideration of loop corrections will be the 
crucial link in this connection. 

In order to provide the reader with some motivation of what follows in the 
course of our report, we shall in this introduction recall the state of affairs in 
topological soliton models, outline the problems one has to face in trying to 
go beyond this status quo, and shed some light on the relationship between 
(baryon) ChPT and the soliton model calculation we are going to consider. 
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1.2 The Skyrme model 



The Skyrme model ||, invented about 35 years ago, yet reconsidered only 
within the timeframe mentioned in the beginning, emerges from the realiza- 
tion that nonlinear effective meson theories allow for topologically nontrivial 
static solutions (solitons) whose conserved topological index is then tenta- 
tively interpreted as the baryon number. Spin and isospin of the soliton are 
supplied by quantizing the zero modes of the problem. The reasons which 
left this model dormant for about 20 years after its inception (except for it 
being somehow orthogonal to the mainstream developments in medium en- 
ergy physics) were essentially twofold: Firstly, in the absence of a suitable 
expansion parameter there was no compelling reason why zero modes should 
be preferred over other fluctuations in the soliton's quantization. Secondly, 
there was no argument to constrain the soliton to have half integer spin. The 
solution to both these problems is of course well known: The difficulty with 
half integer spin could be removed upon consideration of the Wess-Zumino- 
Witten anomaly, whereas the work of t'Hooft 0, |10| and Witten [|l]J on the 



expansion of QCD in powers of the inverse number of colours 1/Nc, besides 
making the concept of solitonic baryons more palatable as a low energy ap- 
proximation to QCD, showed that the pion-baryon coupling is suppressed in 
1/Nc (pions being identified with non-zero mode fluctuations of the soliton), 
and in this guise saved the zero mode quantization as a first approximation 
for baryon properties. 

Since we will be concerned with corrections to this approximation, it is 



at this point worthwhile to recall some results given in reference There 



—1/2 

it was shown that the QCD coupling constant scales like N c if one as- 
sumes a smooth large Nc limit for the 1-loop gluon self energy graph which 
would give rise to the worst divergences in Nc- Essentially from that result 

1 /2 

it follows that the pion baryon coupling must be of 0{N C ), whereas the 
baryon mass is of O(Nc). Provided one accepts the conclusions of [ITT]], one 



can calculate many observables to leading O(Nc), and since the picture of 
meson physics predicting baryon physics seemed rather appealing, the cited 
references consequently triggered a plethora of papers on the Skyrme model 



and its generalizations, examples of which are |12|]-[£6|]. The leading contri- 
butions in 1/Nc to static baryon properties, formfactors and pion nucleon 
scattering were calculated. What marred these calculations and what made 
the 'model' character of the Skyrme model was the difficulty to establish a 
firmly posed lagrangian from which to start, since the 1/Nc expansion gives 
no argument as to the general form of L. This meant that terms in C were 
mostly chosen for computational convenience and carried adjustable param- 
eters. The resulting possibility to 'tailor' the lagrangian in order to obtain a 
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best value for some quantity under consideration makes a clear judgement of 
the merits and defects of the Skyrmion approach to baryon physics difficult, 
because most calculations were limited to one or a few baryon observables 
for which the model used could be fine tuned. Since we want to provide the 
reader with a motivation for our project, we shall now investigate in some 
detail the status of several quantities within the Skyrme model. 

In the following, we present a list of quantities which have been calcu- 
lated in SU(2) Skyrme-type models. Although even this fairly long list is 
not exhaustive, it may serve as a yardstick against which the quality of the 
model can reliably be measured, and which clearly illustrates the problems 



mentioned above. Therefore, in table 1.1 we present entries obtained from 



the standard Skyrme model with one free parameter e, in one case adjusted 
to give the best possible overall picture of baryon properties and in a second 
case chosen in order to cure already at leading O(Nc) what was perceived to 
be a major shortcoming, namely the too large soliton mass. It should be ev- 
ident, that in the latter case, optimization of one value can only be achieved 
at the expense of all others. Eventually, we have included in table |1.1| ranges 
of the results for the different quantities as can be found in the skyrmion 
literature. 



Table 1.1: Tree contribution to various quantities for 
the standard Skyrme model with various values of the 
Skyrme parameter e. Contained is one column that 
lists typical ranges for these quantities as found in the 
skyrmion literature. 





e = 4.25 


e = 7.3 


Skyrmion 
Literature 


exp. 


M[MeV] 


1648 


941 


1300-2100 


939 


a [MeV] 


42 


11 


30-60 


45 ± 7 


<r 2 > s [fin 2 ] 


.98 


.56 


~ 1 


1.6 ± .3 


9a 


.90 


.32 


.8-1.1 


1.26 


<r 2 > A [fm 2 ] 


.09 


.06 


~ .1 


42 + ' i8 

-.08 


<r 2 >E [fm 2 ] 


.61 


.47 


.5-.8 


.59 




1.6 


.47 


1.5-2.0 


2.35 


V? 


.18 


.29 


.15-.25 


.44 


<r 2 > v M [fm 2 ] 


.70 


.40 


.5-.7 


.73 


a [10~ 4 /m 3 ] 


17.0 


2.3 


14-22 


9.5 ± 5 



Some examples from this list are 
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the nucleon mass as a prominent failure of the Skyrme model, since 
it comes out too large by a factor of almost 2 regardless of the detailed 
lagrangian used. An early effort to reproduce the experimental value 



in leading 0(N C ) |T2| required the reduction of f w to one half of its 
experimental value. 

the cr term, where the picture is less clear. Depending on the model 
used, it is over or underestimated. However, this quantity, extracted 
from the isospin even part of the 7rN scattering amplitude has also 
seen a considerable reduction of its experimental value according to 



the latest cr-term update |2"7 ]. The comparably large error bar is partly 
due to inconsistencies in the available data and partly due to some 
model dependence in the data analysis. 

the axial vector constant which is a weak point of Skyrme type 
models, too. It is an excellent example for the statement that adjusting 
the model to reproduce one quantity exactly is very much of a mixed 
blessing with respect to all others. 

the electromagnetic formfactors, where it has long been recognized 
that magnetic moments generally come out too small, most notably the 
isoscalar one, which is wrong by a factor of 2, but also the isovector mo- 
ment for which our parameter gives 1.6 nuclear magnetons in contrast 
to the experimental finding of 2.35. Interestingly, the defective /j, has 
been impossible to repair even with some weird choice of parameters. 

the static electric polarizability, which is generally too large if one 
tries to obtain an overall good fit to the data for other properties as well. 
However, there has been some confusion in the literature concerning the 
exact status of this quantity in soliton models which we will have to 
discuss in more depth later on (section 3.5). 



From the preceeding discussion, two points should be obvious: There is 
a definite need for improvements in the numerical values of several static 
quantities if the soliton picture is to proceed from a qualitative to a more 
quantitative description of the baryon, and such improvements can only come 
from the consideration of higher orders in 1/Nc- Because the next to leading 
order contribution in 1/Nq is carried by pionic fluctuations, this involves 
calculation of pion loops, with all associated problems of regularization and 
renormalization of what is, strictly speaking, a nonrenormalizable theory. 
The second point is that it would be highly desirable to reduce the amount 
of arbitrariness in the model, preferably by establishing a tighter connection 
to meson physics. 
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On the other hand, in view of the obvious problems of Skyrme-type mod- 
els, one might ask why a deeper investigation of the soliton picture is at all 
worthwhile, since HBChPT seems to provide a systematic way to calcu- 
late baryon properties. This point of view, with which we definitely disagree, 
nevertheless deserves a closer look. 



1.3 Baryon ChPT vs. soliton picture: a com- 
parison 

As mentioned above, ChPT, for mesons as well as baryons, is an expan- 
sion in powers of external momenta and quark masses (plus inverse baryon 
mass in case of HBChPT) without any reference to Nq as an expansion 
parameter. Nevertheless, the Nc counting rules outlined above apply gener- 
ally, not just in solitonic theories, and one could look into their consequences 
for theories coupling external baryons to pions . This in fact has been ac- 
complished |p8 |, [p9| . The picture that emerges is the following: Suppose one 



starts with a meson-baryon lagrangian of the form £ = £ M + £mb where 
the subscripts denote the purely mesonic part and the meson baryon inter- 
actions, and wherein the baryon has been introduced as an external field 
and provided £m supports a stable soliton, then the leading 0(N C ) of this 
solitonic theory corresponds to the set of all multi loop graphs without a 
closed meson loop generated from £ and evaluated in the limit of a static 
baryon. (As a shorthand, and in order to maintain the distinction between 
solitonic theory and HBChPT, we shall furtheron call ChPT loops the ones 
generated by £.) Although, from refs |28| ||29|| , this connection is not strictly 
proven, it nevertheless seems very plausible. With the same comment, one 
might go one step further and state that the 1-loop level of the soliton model 
corresponds to the set of multi-loop graphs deduced from £ containing one 
closed pion loop and evaluated in heavy baryon approximation. 

The calculational scheme provided by (baryon) ChPT now ties the number 
of loops to which one carries the calculation to the specific form of £: The 
more ChPT loops, the more complicated becomes £. Stated differently, once 
one has chosen £, one has a limit on the number of ChPT loops to which the 
calculation might be consistently extended. For an infinite number of loops, 
£ would have to contain an infinite number of terms. 

Two consequences follow from these remarks: Provided the picture about 



the connection of HBChPT and soliton models outlined in |28l] |29[ is correct, 



the leading Nc part of the ChPT 1-loop contribution must be contained in 
the tree level of a soliton model. We shall demonstrate this in the course 
of our report on several occasions. The other consequence is that a soliton 



8 



model must be accompanied by an additional assumption not inherent in 
ChPT to justify working with a lagrangian that consists of a finite number 
of terms despite summing an infinite number of HBChPT loops already at 
leading 0{N C ). 

The practicality of a project involving the calculation of loop corrections 
to baryon observables in soliton models therefore hinges upon the answers to 
the following questions: 

• Can we justify the choice of a particular chiral lagrangian used to con- 
struct the soliton, which of necessity must be truncated to some finite 
chiral order, by arguments other than small-derivative-ones? 

• Once this general form of lagrangian then has been established, can we 
make contact to meson phenomenology upon using experimental LECs 
within the theory, and have good values for all baryon properties? 

In the course of this report we shall try to answer these questions. In or- 
der to maintain the tightest possible connection to the meson sector, we will 
show that it is necessary to use a lagrangian that is different from the ones 
found in the skyrmion literature. Before entering the calculation of loop cor- 
rections, we have to recalculate the entire set of tree level quantities from this 
lagrangian, mainly because all tree quantities should be calculated using the 
same lagrangian and the same parameter combination, and because we found 
several earlier tree level calculations to be in need for improvements, specif- 
ically those for the magnetic polarizability and the electromagnetic ratio of 
the A photodecay amplitudes. Additionally, we calculated the neutron pro- 
ton split of the electric polarizability, which had not been derived in soliton 
models before. 

For these reasons and also in order to clarify the relationship between 1- 
loop calculations in the meson and the soliton sector, large parts of the paper 
are of review character. 

This report is organized as follows: Chapter 2 provides the general frame- 
work for the calculations, namely a brief review of ChPT (section 2.1), the 
quantization of the soliton and its renormalization in the presence of those 
external fields related to the baryon quantities of interest (section 2.1) as 
well as the set up of the lagrangian (section 2.2). Chapter 3 then presents 
the investigations for specific baryon observables, in detail mass (section 3.1), 
scalar properties (section 3.2), axial properties (section 3.3), electromagnetic 
formfactors (section 3.4) and polarizabilities (section 3.5) and electromag- 
netic properties of the A resonance (section 3.6). Results are summarized 
and discussed in chapter 4. In this chapter we also draw some conclusions 
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concerning several problems which occured during the process of the 
lations to finally return to the questions posed in the beginning. 



10 



Chapter 2 
General framework 



As pointed out in the introduction, effective field theories are designed to 
describe more fundamental theories (which may well be again effective) at 
low energies where only part of the degrees of freedom are important. In the 
standard model these effective degrees of freedom are assumed to be mesons 
constrained by the requirement of global chiral symmetry which is sponta- 
neously broken into a non chirally invariant vacuum state. Restriction to the 
lowest-lying meson states which may be interpreted as Goldstone modes of 
the broken symmetry leads to the lagrangian 

C eff (U) = £ (2) + £ (4) + . . . £ {N) + ... (2.1) 

where the ChO N counts the number of gradients contained in each term 

A (A ° of 

£W = £ 4 W £ fO (2.2) 

i 

which comes with an LEC to be fixed from experimental information. For 
convenience U is here an SU (2) matrix 

U = e iT77/f * e SU(2) (2.3) 

which depends nonlinearly on the pion field. The restriction to pionic degrees 
of freedom is of course meaningful only below energies where the next higher 
resonances become important. For example, starting from a theory which 
includes explicit vector meson degrees of freedom we could replace the vector 
through pseudoscalar degrees of freedom by virtue of a gradient expansion in 
powers of the inverse vector meson mass my. This would formally generate 
no terms not already present in (|2.1| ), yet make sense numerically only if the 
gradients involved were smaller than my. 

A loop expansion of a theory like ( |2.1| ) necessarily brings about the need for 
renormalization. Provided the regularization procedure used respects chiral 
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symmetry, then ( |2.1| ) is renormalizable in the sense that no counterterm can 
be produced which is not present already. However, this statement is rather 
academical because for practical reasons it is impossible to treat £ e // as a 
whole, instead it must be truncated to some finite number of terms preferably 
at some finite chiral order. The question is of course whether this truncation 
can be justified. 

For the answer, we have to distinguish the vacuum sector from the soliton 
sector. Generally, Weinberg realised that matrix elements of ( j2.1[ ) behave 
like some power p of external momenta carried by the source fields. In a 
graph containing vertices with N derivatives this power is related to the 
number of loops 

p = 2n L + 2 + Y,n N (N -2). (2.4) 

N 

One loop graphs from the N£a model always come along with a power 
p = 4 irrespective of the number of their vertices and the divergencies can 
always be absorbed into the tree graph coefficients of . Proceeding to 
higher chiral orders this situation changes fundamentally. For example al- 
ready one loop graphs from produce all powers p = 4 + 2n^ depending 
on the number of vertices n±, and consequently their renormalization affects 
the coefficients of tree graphs to all chiral orders in (|2.1| ) . 

The characteristic feature of the meson sector now is that external mo- 
menta can be made small by construction i.e. by designing the experiment 
such that this requirement is fulfilled. The truncation of £ in fl2.ip is then 
justified, with ChPT being the result. The crucial difference in the soliton 
sector is that the soliton itself constitutes "external" fields, which cannot be 
made weak by assumption. In fact, gradients of the soliton profiles are typi- 
cally of the order of 700 MeV which is dangerously close to the scale of my. 
Therefore, in the soliton sector, we can neither disregard the counterterms 
of any p generated by a 1-loop graph nor can we, at first glance, get rid of 
multi-loop graphs. 



For the latter problem there is a way out. Following Witten ||11|| , higher 
loops are suppressed by additional powers of the inverse number of colours 
Nc- Thus, for a first correction to the leading tree approximation in the soli- 
ton sector, we consider the 1-loop contribution to be satisfactory. In contrast 
to this, the problem of counterterms to all chiral orders may only be solved 
by an ad hoc assumption, namely we must assume that the renormalized 
LECs of higher chiral orders are small. Whether this assumption is justified 
will be investigated during the course of this report. 
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2.1 Quantum corrections versus 1-loop cor- 
rections in the soliton sector 



As discussed above, in the meson sector, the calculation of quantum correc- 
tions, next-to- leading order in N c , corresponds to evaluating the set of 1-loop 
graphs generated by the chiral lagrangian. This is different in the soliton sec- 
tor. There, the existence of a nontrivial static solution, the hedgehog 

U = e lTfF ^ (2.5) 

means that the classical "vacuum" field configuration is no longer invariant 
with respect to the rotational and translational symmetries of the lagrangian. 
Instead, one has a set of degenerate "vacua" connected to each other by these 
global symmetry transformations. 

A fluctuation which parametrizes such a transformation doesn't remain 
small in the course of its time evolution and therefore has to be treated to all 
orders. Technically, this amounts to introducing collective coordinates [[ItJ . 
It was one of the central assumptions in the early paper on baryon properties 
by Adkins, Nappi and Witten |T2|] as well as in Skyrme's original work || 



that the quantization of these collective modes already supplies a reasonable 
approximation. In this sense we should distinguish in the soliton sector the 
quantum corrections due to collective degrees of freedom, which have been 
taken into account from the earliest attempts on, from quantum corrections 
caused by 1-loop diagrams in which interest emerged much later |^T[, |JT 



2.1.1 Collective coordinates 

As a prelude to the calculation of 1-loop corrections, we shall briefly reca- 
pitulate the collective coordinate method used in the soliton sector |T2]] . We 
introduce the position X(t) of the soliton centre and an SU(2) matrix A(t) 
which parametrizes the soliton's isorotation through three time dependent 
Euler angles 

U(x, t) = A(t)U (x - X(t))A^(t). (2.6) 

Due to the peculiar structure of the hedgehog solution ( |2.5|) a rotation in 
isospace is equivalent to a rotation in coordinate space. With the angular 
velocities AM = \tVL r one obtains 

U = A tym* U ] - XVf/o) ^ (2.7) 
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and the lagrangian 

L = -M + ^M X 2 + ^&(n R ) 2 (2.8) 

in terms of the soliton mass M and the moment of inertia to be specified 
later. For the derivation of fl2.8|) slow motion of the soliton has been assumed. 
With the angular and linear momenta 

R a = 



L a = D ab R b , D ab = ^(r a Ar b A^ (2.9) 
P 91 

1 a 



dX a 



( ( ) denotes the trace in isospace ) following the canonical quantization 
procedure 

[R a , Rb] = —itabcRc, [L a , L b ] = ie abc L c , [P a ,X b ] = —i5 ab (2.10) 
the hamiltonian is obtained as 

H = M + + — (2.11) 

2M 26 V ; 

and for the construction of eigenstates with definite spin (S, S3) and isospin 
(T, T 3 ) quantum numbers 



< x,A\^ s s 3 (P,t) >= ^(-) r -^ ft (Ay( M - a ) (2. 



12) 



SU(2) D-functions may be employed. The calculation of matrix elements 
with the states ( |2.12| ) for the various currents of interest will be explained in 
the corresponding sections. 



2.1.2 Continuum contributions: The phaseshift for- 
mula 

In order to quantize the continuum modes of the chiral lagrangian we consider 
the path integral |37j 



W = N d[U]e- iZ[u] , Z[U] = f d 4 xC (2.13) 
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where Z is the effective action and N a normalization factor. Fluctuations 
are now introduced through the following ansatz |4| 

U = Aie iT ^ x - X)lf ^A\ £ = sju . (2.14) 

A and X denote the collective coordinates introduced in the previous sec- 
tion. The next step then is to expand £ to quadratic order in the fluc- 
tuations. In contrast to these fluctuations which come only with a factor 
I/Ztt = 0(N C 1 ^ 2 ), time derivatives on collective coordinates count as 0(N~ l ) 
and may be neglected because first of all we are interested in the leading loop 
corrections only. Using this adiabatic approximation the path integral reads 

W = N J d[A)d[X)d[ri)e- i I d4x[LiUo ' A ' X)+ ^ Abr ' b] . (2.15) 

For the decomposition of the measure d[U] = d[A]d[X]d[rj\ , orthogonality 
of the fluctuations on the collective modes has to be presumed in order to 
guarantee independence of the integration variables and also to avoid dou- 
ble counting. For that purpose, the zero modes have to be excluded from 
the space of allowed proper fluctuations. Formally, integrating over 77, the 
generating functional 

Z = J d\C(U , A,X) + ^J d 4 x(£nQ) (2.16) 

contains a trace log of the operator governing the time evolution of the fluc- 
tuations . The task now is to evaluate this trace log and, in particular, to 
isolate the divergencies residing in it. 

Here, one encounters the principal technical difficulty of a soliton sector 
calculation: The presence of the static soliton implies that axial, vector and 
(pseudo) scalar sources are contained in the equations of motion (e.o.m) for 
the fluctuations even in the absence of true external fields. These sources 
lead to a metric appearing in the kinetic part of the e.o.m 

n ab = -d t n 2 ab d t -h 2 ab . (2.17) 

Although the fact that the norm n 2 ab is time independent simplifies matters 
somewhat, so far, we are not aware that anybody has been able to perform 
a heat kernel expansion for this kind of operator, which is the method used 
in Q to find the residues of the poles in the trace log of fl. 

Consequently, one has to resort to a different procedure. The alternative is 
the numerical determination of the divergencies from an exact diagonalization 
of the e.o.m for the fluctuations, which according to their time dependence 
~ e~ wt may be written as 

h 2 ab Vb = u) 2 n 2 ab r] b . (2.18) 
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For convenience fluctuations normalized with respect to a flat metric may be 
introduced 



hlbVb = u% , ~h 2 ab = njh 2 cd nj , fj a = n ab r] b , (2.19) 

the spectrum, or, more precisely, the eigenvalues and phaseshifts, are not af- 
fected by this transformation. If h 2 is time independent and does not contain 
time derivatives, one has 

- J d 4 x(£ntt) = - - J d\ (h) , (2.20) 
where T limits the time integral. For a one-dimensional potential scatter- 



ing problem, it can be shown |38[ that the continuum part of the scattering 
operator's trace may be accounted for through a phaseshift integral. Using 
the hedgehog ansatz in the adiabatic approximation the e.o.m for the fluc- 



tuations ( |2.18| ) may be decomposed into partial waves |39j for which the 
one-dimensional phase shift formula applies. The total trace of the scatter- 
ing operator is then obtained by replacing 5 by the sum over phaseshifts of 
all angular momenta 

l r ,r . , l r°° . d 



- [ d\ (h — ho) = — / du uj ——5{uj 
2 J N u/ 2tt Jm n du K 



1 f°° i d 

J dpy/p* + ml T 6(p). (2.21) 



2tt Jo w dp 
Here, h\ = —A + m 2 corresponds to the free Klein-Gordon equation in the 



absence of the soliton fields. In the last step the u> integration, u> = \jp 2 + m 2 , 
has been converted to a momentum integration. From (2.21) it is clear that 



divergencies are related to the high momentum behaviour of the phaseshifts 

S( p ) ^ a p 3 + a lP +^ + 0(p- 3 ) . (2.22) 

V 

The explicitly denoted terms give rise to at least logarithmically divergent 
expressions. The strategy is now to subtract the worrisome terms in the phase 
shift integral and add them separately. Using dimensional regularization 



dd V ,_2 , _2M/2 _4 f, , 1 / 1 | p™l 

\2n* [ 2 + p? 

d "*"?? , o o,_i/o . o I , 1 TTl 2 .\ 



(2vr) d - 1 ^ 7Tj \ 32vr 2V /i 




,2 



dd ~ lp -(p 2 + ml)-^ = -8(\ + ^(1 + ItJ%)) , (2.23) 
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which involves a scale to render the arguments in the logarithms dimen- 
sionless, the poles as d — > 4 which reside in 

A =£(^4-^( 1 ) +fc ^' +1 < 2 ' 24 > 

may be isolated 

]-Jd\(h- h ) = E cas {fi) + A(37rm^a - Anm 2 ^ + 87ra 2 ) . (2.25) 
The Casimir energy 

r°° dp 

JO y^2 _|_ m 2 

+ ^l^.^f^e^ ( , 26) 



1 
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= - a p 3 - dip) + a 2 ] - m n 5(0) 



is a finite but scale-dependent expression. For the derivation of eqs.( [2.25" 



2.26| ) an integration by parts has been performed and m 2 has been added 



to the denominator of the a 2 term in order to avoid infrared problems. The 
coefficients ao,ai and a 2 contain in general all ChOs as is noticed from the 



fact that h = Jho + w with some potential w. Therefore the trace ( [2.25] ) 



contains not only all ChOs via the soliton's stability condition but also an 
infinite numer of explicit gradients. The part proportional to the divergent A 
may formally be expanded 

3vrmtao - ^m\a x + 8vra 2 = £ i\ N) j d\ C\ N) , N > 4 (2.27) 

i,N 

into the complete set of terms contained in the chiral lagrangian ( J2.1| ). 
The coefficient 7- for a definite term of ChO N is in general complicated 
and scale dependent through all the renormalized LECs of ChO smaller 
than N. Again, it is only the 1-loop contribution to the N£a model which 
corresponds to the lowest no n- vanishing ChO 4 in ( |2.27| ), where the situation 
becomes simple. For that particular case (a^ = 0) and a 2 are ana- 

lytically known M p4 and the coefficients 7$ are simple numerical factors 



(section 2.2). 

According to Q2.27 ) the divergencies in ( |2.25 ) may finally be absorbed 



into a redefinition of the lagrangians LECs. The total soliton mass (tree + 
1-loop) is obtained as 



M (£\ N)r ) + E C M , (2-28) 
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where E cas represents the Casimir energy (|2.26|) and 

M^f )T ) = -J2 £f )r [ d 3 r £™ (2.29) 

i,N J 

is the tree contribution to the soliton mass calculated with the renormalized 
LECs 

^ = 4*° - T^A . (2.30) 

It is this relation which for N = 4 establishes the close connection to ChPT 
in the meson sector with the conventions adopted by Gasser and Leutwyler. 
With A defined in ([2.24 ) the renormalized LECs £^ r (|2.30|) just coincide 



with those defined in M and auspiciously we may take over their analysis 
(section 2.2). 

According to ( |2.30|) the renormalized LECs for the higher ChOs, in par- 
ticular their scale dependence, are fixed in principle, however for practical 
purposes we do not know the 7,- for the individual higher ChO terms. Even 
if we knew the value of the corresponding renormalized LEC at some scale 
we could not, according to ( |2.30|) determine its scale dependence without 
knowledge of the Ji ■ For a detailed examination of the scale dependence of 
the soliton's mass see section 3.1.2. 



As mentioned, the phaseshift formula (|2.21|) takes care of the continuum 



only; if there exist any bound states at energies u c the zero point energy 
\ He w c has to be added to the Casimir energy. In the absence of external 
fields, the spectrum of h 2 contains no true bound states, instead there exist 
zero modes due to the rotational and translational symmetries discussed in 
section 2.1 . Because of u c = the zero modes do not explicitly contribute. 
However, in the presence of external fields, which violate rotational or trans- 
lational symmetries, zero modes are shifted to finite energies. Their treatment 
is described in section 2.1.4. 



2.1.3 External fields 

The above formulae allow to calculate the static soliton's energy including 
quantum corrections due to pion loops. However, the aim of this investiga- 
tion is to calculate other baryon properties on the same footing, too. Such 
quantities may be defined as the linear change of the baryon's energy in the 
presence of an external stimulus. 

To be specific, we study the coupling of an external field j with strength 
e to the lagrangian fl2.1|) 

C(e) = C{i h U) - ej ■ J{£ h U) , (2.31) 
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where J(ii, U) denotes the corresponding current density. The external field 
has to be chosen suitably so as to give the desired quantity, e.g. for the cal- 
culation of the magnetic moment, J represents the electromagnetic current 
density and j ■ J corresponds to the magnetic moment density. Matrix ele- 
ments of j ■ J are then obtained as a derivative of the soliton mass (tree + 1 
loop) in the presence of the external field with respect to its strength 

(N\fSrjJ\N)^ 9M{£) 



9S e=0 

M(e) = M (e,£ r l )+E cas (e,fi) ■ (2.32) 

Here, \N) represents the one baryon state ( |2.12|) . 

Considering static baryon properties, the external field can be specialised 
to a simple, in many cases space-time-independent form, and can generally 
be assumed to be weak. Consequently, we shall proceed as follows: 

(i) The external field is chosen such that the quantity of interest is obtained 
from ( p.32 ) in tree approximation. 



(ii) The static soliton is computed from ( |2.31| ) for some small e. It turns 
out that for all quantities considered here, the hedgehog still solves the 
classical e.o.m in the presence of the external field. 

(iii) The operators h 2 (e) and n 2 (e) are obtained by expanding C(e) to 
quadratic order in the fluctuations. The stability condition for the 
hedgehog (ii) ensures that the term linear in the fluctuations is ab- 
sent. 

(iv) The e.o.m for the fluctuations h 2 (e)r] = uj 2 (e)n 2 (e)rj(e) in the pres- 
ence of the external field is solved for the scattering phaseshift up to 
sufficiently large linear and angular momenta, respectively. 

(v) The Casimir energy E cas (fi,e) is computed according to (|2.26| ) and the 
soliton mass in tree + 1-loop is obtained in the presence of the external 
field. 

(vi) Repeating the whole procedure for several small values of e , the deriva- 
tive dM(e)/de\ £=0 may finally be computed, which equals the desired 
quantity including quantum corrections. 

We again stress that in this method the integral has to be restricted to the 
continuum fluctuations orthogonal on the zero modes. Whereas this orthog- 
onality holds automatically in the absence of external fields, it generally has 
to be imposed in their presence (c.f. next section) . 
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A second caveat concerns eq. (|2.20 ) . If the external field creates a time- 



derivative interaction as is the case for the electric polarizabilities (section 
3.5) 

Q, = -d t n 2 d t - 2iewd t - h 2 (2.33) 
where w (and w = n~ 1 wn~ 1 ) is time independent, we have 

i r .a ,„ „, T 



% - J d A x(lnVL) = i - J dp J d\ (£n(p 2 - 2ep w - h 2 )) 

J d\ (^h 2 + e 2 w 2 ) (2.34) 



T 
2 



i.e. there is no term linear in e and the leading contribution is of 0(e 2 ). 

The specific choice of external field, as well as other details, will be given 
separately for each quantity in question within the corresponding subsection 
of the next chapter. 

2.1.4 Treatment of zero modes in the presence of ex- 
ternal fields 

Due to translational and rotational invariance the system of e.o.m for the 
adiabatic fluctuations possesses 6 zero mode solution z c a with oj c = 

/&z 6 c = 0, c=1...6 (2.35) 

which enter the phaseshift formula for the Casimir energy only indirectly via 
Levinson's theorem 5(0) = 6n. Because of u c = , there is no direct bound 
state contribution as e.g. in the case of the kink's breathing mode [30|. 



In the presence of external fields, the situation changes, since rotational 
and/or translational symmetries are in general violated. This is always the 
case if the external field coupled to the lagrangian is not rotationally or 
translationally invariant. If such a field is switched on weakly with strength 
e, the zero modes of the corresponding broken symmetry are shifted to finite 
energies u c ~ yfe. This is understood immediately, because the external field 
generates an additional term in the e.o.m for the fluctuations h 2 (e) = h 2 + ew 
which in perturbation theory leads to u 2 = e J dh" z^v a bZ^ if the matrixele- 
ment does not vanish because of rotational or translational invariance. For the 
case of the axial vector coupling constant gA and the corresponding external 
axial field which violates the former but not the latter symmetry (c.f. section 
3.3 ) we calculated the energies of the rotational zero modes numerically for 



different field strengths e. The result, plotted in fig. 2.1 illustrates the be- 



haviour u c ~ ^fe for the zero modes which, considered as true bound states 
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with contribution ^ J2 C w c to the Casimir energy would lead to a desaster: The 
corresponding baryon property dM(e)/de\ £=0 would aquire infinite quantum 
corrections. 



Figure 2.1: Energy of the rotational zero 
modes as function of the strength e of a con- 
stant external axial field. Because the external 
field violates rotational invariance, the zero 
modes are shifted to finite energies propor- 
tional to the square root of the external field 
strength e. 



The solution to this problem is obtained by removing the zero modes, 
which are not confined by a restoring force, from the space of allowed small 
amplitude fluctuations. In fact, their contribution has already been taken into 
account in tree approximation by the introduction of collective coordinates. 
To remove the redundant variables in order to avoid double counting we must 
impose constraints JHJ (see also section 3.5 for details) 

& z>l b Vb = 0, !d\ z c a n 2 ab z d b = 5 cd (2.36) 



which have to be added with multipliers A c to the lagrangian. Only then 
the Legendre transformation to the hamiltonian will be well behaved. With- 
out external fields, the constraints ( [2.36D are automatically fulfilled and the 



calculation of the soliton's Casimir energy remains untouched. 
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However, in the presence of symmetry breaking external fields, the con- 
straints become active and contribute to the e.o.m 



h 2 ab (e)r] b (e) - X c (e)n 2 ab z c b = u 2 n 2 ab (e)rj b (e) . (2.37) 

The dependence on e is indicated wherever necessary; for the constraint itself 
this dependence is irrelevant, since 0(e) is sufficient to calculate the baryon 
property dM(e)/de\ £=0 of interest and A c is already of 0(e) : 

X C = J d\ z c a [hl h {e) - u 2 n 2 ab (e)}r, b (e) ~ e (2.38) 

In fact, we are going to demonstrate that apart from removing the zero mode 
bound states the effect of the constraints on the scattering states is 0(e 2 ) 
and may therefore safely be discarded. 



For this purpose, we first solve the e.o.m (|2.37| ) without constraints 

h 2 ab (e)f) b (e) = uj 2 n 2 ab (s)fj b (e) (2.39) 

and then treat the constraints as perturbation calculating the DWBA matrix 
element which is related to the phaseshift change 

A5(p) =p{jd\ fj a n 2 ab (e)zt] {/ d\ z c e [h 2 ed (e) - co 2 n 2 ed (e)}Vd} ■ (2.40) 

Both integrals vanish linearly with e and the phaseshift change is clearly 
0(e 2 ). 

Thus the situation becomes quite advantageous in this respect: For the 
phaseshift calculation the constraints may be ignored and in the Casimir 
energy they just remove the unwanted contribution \ Y JC ° J c- The phaseshift 
formula ( 2.25 , 2.26| ) as it stands remains correct in the presence of external 



fields which violate rotational and/or translational symmetry. 



2.2 Chiral lagrangian in the soliton sector 

We now turn to the discussion of the lagrangian to be used in the soliton 
sector. Since for consistency reasons we have to solve the static e.o.m of the 
same lagrangian that generates the loop graphs, it is clear that such an object 
has to contain at least ChO 4 terms - an N£cr model doesn't support a stable 
soliton. The question then is, what chiral order, if any, would be sufficient? 

Recalling the remarks of section 2.1 we have to face the facts that coun- 
terterms will contain all chiral orders and, even more disturbingly, that an 
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expansion in powers of external momenta must fail because of the size of 
gradients on the soliton field. 

Consequently, it would be meaningless to simply count gradients in 
Weinberg-Gasser-Leutwyler fashion. Instead, either the renormalized LECs 
beyond some chiral order must be small or a cancellation between different 
terms must occur in order to assure negligible contributions. 

The smallness of a renormalized LEC is of course to some extent depen- 
dent on the choice of scale. But this choice is restricted in several ways: A 
pragmatic restriction stems from the requirement of a small symmetric 4th 
order term in £ so as not to destroy the soliton. Another, more fundamental 
point is that the renormalized LECs must be adjusted to be predominantly 
of 0(Nq) since otherwise the 1/Nc expansion must fail. The lagrangian up 
to ChO 4 given by Gasser and Leutwyler M reads 

C GL = Jl <a ^>+^-<U + U^> 
4 ^ 4 



+ j < > 2 + j < a^a u 



+t\ < F^ V UF^ R U ] > 

+*■! < F^ u D»UD v U ] + F^D fl WD u U > , (2.41) 



where 



U^D^U, = d ll + i[v ll ,-] -i{a^-} 



F nv = 9f,(vu ± a v ) - d v {v^ ± a M ) - i[(v^ ± a^), {v v ± a v )\ 

and f M (a M ) are external vector (axialvector) fields. Renormalized LECs and 
the respective 7« are given in table |2l| Using these numbers ( [2.41|) allows for 



a soliton solution; however, with the value of l r 2 corresponding to a Skyrme 
parameter e = 7.24 this soliton is unphysically small, compare table pi 
Therefore one has to conclude that something is missing. To quantify 'some- 
thing' and to get a better understanding of higher ChO contributions at least 
at tree level, we investigate in the following a pionic N£a model coupled to 
several other mesonic degrees of freedom, namely scalars, isovector vector 
and isoscalar vector mesons. 

Neither is this list of resonances exhaustive nor is their inclusion into 
a chiral lagrangian unambiguous. Concerning these objections, at least the 
quality of the vector meson dominance (VMD) assumption on which our 
way to incorporate vector mesons is based can be judged using experimental 
values for ChO 4 LECs reference. 
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2.2.1 Contributions from scalar mesons 



To estimate the effects of scalar mesons in a purely pseudoscalar model, we 
start from an N£cr model in which a scalar a meson has been introduced as 
a dilaton [TJ5[], JIT] fj"2"fl . This model contains two parameters, the glueball 
(SOOMeV") 4 [[|3|] and the glueball mass m G 
\2mMeV. 



condensate Cq 



1200MeV 



leading to a scalar meson mass m a 

In order to obtain a purely pionic lagrangian up to ChO 6, we expand the 
scalar a field in powers of 1/Cg- Only the first two terms in this expansion are 
of relevance, higher ones would affect the result only at ChO > 6. Inserting 
the approximate a field back into C, we arrive at 



C 



r 



f 4 

J 7T 



{a^ - 2ml(l - u)} 



2 

+ 



(2.42) 



mhCG 1 ' 



where we used 



a 



(2) 



a. 



II - 3ml(l - u)} 

it ■ a M , u = - < U + > 



[2A3) 



Numerically, we have calculated the chiral angle F from the original la- 
grangian [19| and inserted into ( p.42j ) (note that for the present purpose 
( [2.42 ) is understood not to contain external fields). The contributions to 
mass, ttN a-term and qa are then compared in table 2.1. Convergence in 



this case is excellent. (The Skyrmeterm, which is not affected from a meson 
exchange was omitted in this comparison.) 

Likewise LECs could be read off and confronted with experimental values; 
however, we cannot expect agreement since we deliberately chose only to 
include the dilaton a meson; therefore we skip this here. 



2.2.2 Contributions from vector mesons 



Proceeding to calculate in analogous manner the effects of vector mesons, 
we use as a starting point the so called minimal model [I3| including 
symmetry breakers |44| and coupled to external sources. This is not the most 
general implementation of uj and g mesons as far as anomalous pieces of C 
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Table 2.1: Comparison of exact and approximate scalar 
meson contributions. Each of the three intermediate 
columns contains the sum of approximate pieces up to 
and including the indicated ChO. 





ChO = 2 


ChO < 4 


ChO < 6 


exact 


Mass [MeV] 


818.5 


691.7 


712.8 


716.9 


a nN [MeV] 


45.8 


28.9 


34.0 


34.4 


9a 


.473 


.384 


.40 


.404 



are concerned jy| [Q , but the differences can be shown [[|6|] to manifest 
themselves only at ChO > 8. 

Expansion of the g and u fields in powers of the inverse vector meson mass 
and subsequent reinsertion into C yields (again up to ChO 6) 

f 2 



£ = £ {a^- 2rn*(l -«)} 



2 

\ < effi" > < - ^ $])r"« > -^b* 



2m 2 v ^ '^ J ' 2 

' < P'VS - ig[v^ &f > -P= < r?V*« > 



Ami " 4m ' 

- ^> M3) + 2(^fB„B» . (2.44) 

Here, the following abbreviations have been used: 

( 3 ) N c g„ 
mi 

? ] = A t ([C/ ) « M ] + [a M ,C/t])e. (2.45) 



r 



_B M is Skyrme's baryon current. The strength d of the standard symmetry 
breaker for vector mesons is related to the 577(3) meson mass differences 
d = {m 2 K , — m 2 g )/8(m 2 K — m 2 ). Moreover v Q means an isoscalar vector field 
and the field strength tensors are defined in the usual way 

fnv = d^f u - d v f lx - i(g)[f^, f v ] 
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where the factor g is in place if / denotes a vector meson while it is absent 
in case of external fields. The coupling to external fields was considered since 
the local approximation can result in nonminimal couplings. 

As an aside, one notes that go as well as uji are induced through the 
external vector field. Such components are important for certain properties 
(e.g. polarizabilities, see section 3.5) and have so far been missed in the 
calculations. 

Without external fields, numbers are again computed by inserting the 



exact chiral angle obtained from the lagrangian according to [IS| for 



g = 2.9, g u = 2.2 into fl2.44f) . The comparison in table |2.2| then shows that 



the ChO 6 contribution due to the g is indeed much smaller (by a factor of 
4) than the one of ChO 4; nevertheless the approximation does not improve - 
the exact value is, including ChO 6 underestimated by the same amount that 
it was overestimated using only ChO 4. For the uj the first approximation 
differs from the exact result by a factor of 2. 



Table 2.2: Comparison of exact and approximate contribu- 
tions to the soliton mass (in MeV) for g and uj mesons. The 
two intermediate columns contains the sum of approximate 
pieces up to and including the indicated ChO. They are fur- 
ther divided so as to display the contributions to symmetric 
(~ mj) and symmetry breaking (~ mf£ n \ n > 0) pieces 
individually. Note that, although the original Lagrangian 
contains terms with n = 1 at most, its approximate may 
comprise terms of arbitrary n. 





ChO < 4 

n 2 


ChO < 6 

m m 2 m 4 


exact 
n 2 


Q 


494.6 


377.6 5.9 -.36 


439.0 5.66 


UJ 




475.3 


235.85 



Since uj can be expressed exactly in terms of the baryon source, it is easy 
to go beyond this approximation. We have 

- , ^3 e iq{r — r') 
^"ZfinV) (2-46) 



for the (only nonzero) static field ujq. If 8^Bq and are of the same mag- 
nitude, the approximation will not converge. This indeed seems to be the 
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case. 

From the above reasoning, it is obvious that the ChO 4 lagrangian ac- 
counts for all effects from scalar mesons and to a large extent for q mesons, 
too. This last statement is also obvious from the comparison of VMD pre- 
dicted versus experimentally known LECs, table [2.3| . (The ji used in this 
table are taken from [|J, LECs from |47|). 



Table 2.3: LECs at scales of \i 
compared to VMD prediction. 



and /i = m Q , respectively, 





li 


= m v ) ■ 10 3 


= m e ) ■ 10 3 


VMD 


i\ 


1/3 


-3.94 ±1.3 


-4.65 ±1.3 




pr 


2/3 


6.35 ±1.5 


4.92 ±1.5 


l/(16s 2 ) =7.3 


pr 


-1/2 


-0.22 ±3.8 


-0.85 ±3.8 




t\ 


2 


9.74 ±5.7 


5.46 ±5.7 




11 


-1/6 


-5.88 ±0.7 


-5.22 ±0.7 


-1/(16# 2 ) = -7.3 


pr 


-1/3 


-14.5 ±1.2 


-13.8 ±1.2 


-l/(8# 2 ) = -14.6 



Figure 2.2: Comparison of energy densities of 
purely pionic models with the three parameter 
sets A, B and C (solid lines) and the exact 
density from the VMD lagrangian with g = 
2.9, g u = 2.2 (dashed line). 
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What is missing in ( |2.41| ) essentially is due to the to meson. How then can 
we account for this? To answer, we compared the energy density of a VMD 
minimal model with g = 2.9 g u = 2.2 fl8| to its local approximate ( fl2.44j ) 
without symmetry breakers and 6th order contributions from g). The chiral 
angle was in this case computed from the approximate lagrangian, and the 
result is plotted in fig. [2.2| . (For a converging local approximation, this com- 
putational difference shouldn't matter; however, since we suspected lack of 
convergence for u, and since we would use a chiral angle from a purely pionic 
model in our calculations of 1-loop corrections, we had to do the comparison 
in this way.) From fig. it is obvious that much better agreement may 
be obtained by calculating the chiral angle from an approximate lagrangian 
with parameters reduced to g u ~ 1.0 and e = 4.5: the missing higher chiral 
orders renormalize the original parameters to effective values. Eventually, a 
computationally convenient model with e = 4.25 and g u = is not worse 
than the original e = 5.8, g w = 2.2 combination. 

These considerations result in the lagrangian 



-ol^o? + ml(u - 1) 



+ -c s 4 (ac u a 



M2 



+ cl(a^)(u-l)+cT(u-l) 



r f 

+ c\ < F^UF^ uR U ] >+ij< F^D"UD U U^ + F*D»U ] D U U > 

1 -v^B^ (2.47) 



0^ - 2n%B^ 



Am 2 , 



At 2 



+ 



ft fte* 



P 

J TT 



m 



2 *-4 
P 

J TT 



m 



4<-3 



r + 



p 

J TT 



±5. 

p 

J TT 



— c 



To — — TT c 6 



ft 8 ^ A^fl 

Note that only the photon can couple nonminimally to the baryon current. 
For comparison purposes, in the following evaluations, we therefore use ( |2.47| ) 
with the parameter combinations mentioned above, which will henceforth be 
referred to as models 

A = (e = 4.25,^ = 0) , 

B = (e = 4.5,&, = 1.0), 

C = (e = 5.8,^ = 2.2) . 

These parameters are used to determine cl and cq as well as the constants 
cf , c{ regulating the contributions of nonminimal terms, whereas c\ and c™, 
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both scalar meson induced, are kept at their experimental values for /x = m e . 
One notes that the constant in front of the nonminimal term produced by 
the u) depends solely on the vector meson mass and is therefore not affected 
by the use of 'effective' coupling constants. Unless otherwise stated, c| is set 
to zero in accordance with experiment. Finally, the finite renormalization of 
/ — > f n , m — > m n is expressed through the formulae 

fl = f + <lt\m\ flml = fm 2 + 2(*£ + Qm' . (2.48) 



To summarise, we have investigated the local approximations on models 
containing explicit scalar and vector meson degrees of freedom in order to 
better understand the role of higher (ChO > 4) chiral orders in the chiral 
lagrangian, which, in contrast to standard ChPT cannot be dismissed out of 
hand in the soliton sector. We concluded that terms of higher chiral orders due 
to scalar mesons and isovector vector mesons are essentially unimportant, but 
that there is no a priori justification to drop those induced by the isoscalar 
vector meson. Since we could not take those terms into account properly, 
we used, as an alternate means to handle this problem, 'effective' coupling 
constants in the terms generated by vector mesons These coupling constants 
differ from the ones in the original vector meson lagrangian; comparison of 
exact and approximate energy densities led to parameter combination B, 
whereas A and C would serve as reference sets, with C using paramters un- 
changed relative to the exact lagrangian and A being particularly convenient 
to work with. 



2.3 One- loop corrections for lagrangian with 
explicit vector and scalar mesons 

Although we use a purely pseudoscalar lagrangian for the calculation of 1- 
loop corrections in this report, we briefly discuss the inclusion of other me- 



son species like vector and scalar mesons. Such models [I5|]-[|19| have 
been quite sucessful in some respects (high-energy behaviour of pion-nucleon 
phaseshifts, for mf actors) . 

The most general lagrangian of this kind would be of the form 

C x = C eff (U)+J2^nt(U,x)+J2^res(x) , (2.49) 

x x 

where x = g,u,a... generically denotes all possible resonance degrees of 
freedom and C e ff is of the form ( |2.1[) . It is clear that such an object would 
be even less manageable than the purely pseudoscalar lagrangian and has to 



29 



be accompanied by some simplifying assumptions. The first one is obviously 
the restriction to low lying resonances, i.e. vector and (two kinds of) scalar 
mesons. The second one uses the fact that upon expanding the resonances 
in terms of pseudoscalars, C x itself is of the form fl2.1| ). Calculating 1-loop 
processes, we have £^ — > £\ N ^ r (fi) as previously, but for N > A, i\ N ^ r may 
be decomposed into a resonant and a direct part 

4 ^ (A 4 ) = 4 res + 4 direct (A*) j (2.50) 

where £fl es depends solely on resonance mass and coupling constant, but 
carries no scale dependence. The assumption then is that there exists a scale 
such that 4&rerf(/ i ) — 0> which is an alternative way to spell the concept of 
vector (and scalar) meson dominance. This postulate is considerably weaker 
than its analogue in the purely pseudoscalar case, namely to assume the 
existence of a scale where i\ N ^ (fi) ~ for all higher chiral order terms which 
cannot be accomodated within our formalism. This constitutes the main 
conceptual advantage of a calculation involving resonances explicitly. The 
total soliton mass in tree + 1-loop ( |2.28|) generalizes to 

M(ri = M + Y,E x c M- (2-51) 

X 

The classical soliton mass M depends on parameters introduced by the addi- 
tional mesons, e.g. g g , . . . . All mesons contribute to the Casimir energy via 
their phaseshift 5 x (p) which is a sum over all channel eigen phases. The num- 
ber of channels may be considerable, e.g. 9 channels for the p-meson, which 
makes an accurate determination of the phaseshifts technically difficult. A 
nice feature is however that in the high momentum region these phaseshifts 
behave well, i.e. in the individual partial wave they tend to zero in contrast 
to those of purely pseudoscalar models. Consequently the determination of 
the asymptotical constants 

5 x {p) p ^ a x p + ^ + 0(p~ 3 ) (2.52) 

should be less critical compared to the pure pseudoscalar case which involves 
an additional term a p 3 (compare eq. ( |2.22| )). With these preparations the 



contributions of the individual meson species x to the Casimir energy may 
be evaluated with the phaseshift formula 

Ztt 



'''' ' p(S x (p)-afp)+a x 2 }-mJ x (0) 



p 2 + m 2 . 

mltf , ml a\ . ml. 
4 /i 2 2 



(2.53) 
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Because the meson mass m x is smallest for pions it is expected that their 
contribution dominates the Casimir energy, the contributions of other mesons 
being suppressed by their larger mass. 

There remains the question concerning the chiral scale /i. For a it qui model 
which in addition includes both sorts of scalar mesons the LECs in ChO 4 
compare well [H71 with those obtained in ChPT at /i ~ m g (table |2.3| ). 
Therefore, such a choice should be quite reasonable. Anyhow the Casimir 
energy will not react very sensitively on small changes in the chiral scale 
because it enters only logarithmically. 

Below, we list merits and problems of using lagrangians with explicit 
vector- and scalar-mesons. 

Advantages: 

• The LECs of higher chiral order terms are fixed form the assumption 
of vector (scalar) meson dominance and not simply neglected as in the 
pseudoscalar case. 

• the high momentum phaseshifts behave well 

• all parameters are fixed in the meson sector, these models possess no 
free parameters. 

Disadvantages: 

• the inclusion of additional mesons is not unique as far as the anomaly 
is concerned 

• next higher resonances (e.g. axial vectors etc.) are treated only approx- 
imately 

• the tree calculations have to be updated (induced components, e.g. 
polarizabilities) 

• a huge coupled channel problem has to be solved with sufficiently high 
precision. 

In the end we decided to use a pseudoscalar lagrangian only. We are con- 
fident that, concerning 1-loop corrections, the results for a lagrangian with 
explicit q and u -mesons would come close to those of our model B. In par- 
ticular the problems with the axial quantities discussed in section 3.3 are not 
expected to be cured by the introduction of vector mesons. 
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Chapter 3 
Baryon properties 



In the following sections, we are going to display the specifics of the way to 
calculate various baryon properties of interest. 

For the presentation, we shall adopt the following pattern: For each quan- 
titiy under investigation, we will first give some basic definitions and exper- 
imental findings. In parallel, we explain the choice of external field so as to 
admit a hedgehog solution. We then show a recalculation of the tree level 



value comprising all relevant terms of the lagrangian ( 2.47 ). In case of the 
electric polarizability, doubts had been voiced regarding the correct way to 
calculate it; we will discuss this issue and conclude that the standard way is 
indeed correct. 

Proceeding to the 1-loop calculation, we derive the e.o.m for the fluctu- 
ations and discuss problems which arise in several instances. They concern 
the proper meaning of vacuum subtraction (a-term, e.m. formfactors) and 
the treatment of apparently non local terms (electric polarizability). 

For a more convenient comparison, numerical results will then be discussed 
in a separate section. As always, there is an exception from the rule; in the 
present case this involves the axial coupling constant, where we give results 
immediately since we are forced to investigate in detail the implications of 
current algebra for the 1/Nc expansion of this quantity. 

The full list of calculated quantities involves 

• the Baryon mass (section 3.1), 

• the ttN o term and scalar radius (section 3.2), 

• the axial coupling constant g& and the corresponding axial radius (sec- 
tion 3.3), 
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• the electromagnetic formfactors (more specifically, the isovector mag- 
netic moment and corresponding radius and the isoscalar electric ra- 
dius) (section 3.4), 

• the electric polarizability of the nucleon (section 3.5). 

• the electromagnetic properties of the 5 isobar (section 3.6). 

With respect to the polarizability, we opined that its neutron proton split 
deserves a derivation, which has therefore been included in section 3.5.3 . 

Although the list of calculated quantities is fairly long, it is by no means 
exhaustive, since we are limited to quantities where the correction is brought 
about by 'adiabatic loops', i.e., loop graphs, in which the pion can be treated 
in adiabatic approximation. A calculation involving nonadiabatic fluctuations 
would also have to deal with (adiabatic) 2-loop graphs which appear at the 
same level in the 1/Nq counting. Such a project is well beyond our present 
abilities. 



3.1 Baryon mass 



3.1.1 Tree approximation 

The classical soliton mass is obtained by inserting the hedgehog ansatz into 
the lagrangian (2.47). For convenience, we will give this quantity in terms of 
the longitudinal and transversal metric appearing in fl2.17j) . 



n lb = h(r)fah + b T (r)(5 ab - r a r b ) 



(3.1) 



which allows to express the mass and most of the other quantities in the 
subsequent sections in a very compact way without reference to the specific 
terms contained in the lagrangian. For our choice, these functions are given 
by 



b L = l + 2ct S --2ct(F- + 2 -^)- 
b T = l + cl(F' 2 + ^)-2cl(F' 2 + 



2cJ(l - c) + c 6 - 



(3.2) 



2s 2 



F ,2 g 2 



-2c*(l 



where abbreviations s = sin(F), c = cos(F) were used. 
Variation of the classical soliton mass 



c) + Or 



f 2 

y 



F% L + 



2s 2 



(3.3) 
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with respect to the chiral angle yields the stability condition 



\(r 2 F'b L )' 



2sc o 
—b T + m^s 

2 s 2 



2c™s(l 



(3.4) 



an ordinary nonlinear differential equation, which has to be solved numeri- 
cally subject to the boundary conditions -F(O) = 7T, F(oo) = which guar- 
antee a solution carrying baryon number B = 1 . Clearly the soliton mass is 
of0(iV c ). 

It was the fact that the classical mass ( ^.3[ ) always came out much too 
large in soliton models with realistic parameters which finally led to the 
investigation of loop corrections [ET]] [3~I]-[[5E[] which in the end turned out to 
considerably lower |3~" 



[p5 | p6[ the numerical values for this quantity 
Similar to the soliton mass, we obtain for the moment of inertia 

e = Hfi J d \ s% T , 



which according to ( 2.11 ) determines the nucleon-A split 
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(3.5) 



(3.6) 



a nonadiabatic quantity of 0(N^ ) related to the angular rotation. As al- 
ready mentioned, we do not report about loop corrections to nonadiabatic 
quantities although they are by no means less important for the simple reason 
that such calculations would become tremendously complicated. 

Another nonadiabatic quantity related to the mass is the neutron proton 
split which vanishes for the lagrangian fl2.47p because there exists no term 
that distinguishes between states of different isospin 3-component. In fact, 
it can be shown that such a terms must not appear before ChO 8 |4l][] . The 
ChO 8 term is isolated by local approximation from the standard symmetry 



breaker for vector mesons and is related to qui mixing [g4j |^8 



L 



tree 



16g s m 2 g ml 
)D 3a 



> 



(6M, 



rip) 



(3.7) 



Here, is again the familiar baryon current. With experimental values for 
g u = 3.3, g e = 2.8, m u = 782MeV, m e = 768MeV and for the qoj mixing 
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parameter \x 



-4.5 ± 0.6) ■ 10~ 3 GeV 2 p| p| quite a reasonable result 



is obtained for the np split 



M, 



II p 



30g e m 2 g ml 
compare table PO . 



2F's 



+ 



s 2 c, 



B° 



F's 



2ttV 2 



(3.8) 



3.1.2 Loop corrections and scale dependence of the 
soliton mass 

For the 1 loop corrections adiabatic fluctuations around the soliton back- 
ground according to the parametrization ( 2.14j) are introduced and conve- 
niently decomposed into longitudinal and transversal components 



r\ = n L r + rj T . 

From the lagrangian expanded to quadratic order in the fluctuations 



(3.9) 



\JSr {if 



4c 



2(c 2 -s 2 ) 2 

+ ,2 W L 



(F' 2 + KWt + 



[vanl b Vb + VaKbVb] 



(3.10) 



the operators n 2 b and h 2 b which determine the e.o.m ( ^.18[ ) may be read off. 
For simplicity, we give only the terms generated by the N£a model explictly. 

Technically, it is convenient to write the fluctuations in terms of vector 
spherical harmonics 



Ll 



-iu) t 



(3.1i; 



which decouple the e.o.m into electric £ — L ± 1 and magnetic i = L modes. 
This procedure is standard and the full differential equations for the radial 
functions /«(r) are given elsewhere J3"5"| 



pHJ . The challenge is now to 
solve these coupled channel equations for the phaseshifts with high enough 
accuracy which implies large phonon spins (L max ~ 100) and large momenta 
(p max — 25m n ) such that the sum 



S(p)=Y,(2L + l)Sl(p) 

Lc 



(3.12) 



35 



converges and the asymptotic behaviour can be reliably extracted. For this 
purpose we use the so called variable phase method |52| which allows to 
solve for all desired phonon spins simultaneously. The result is plotted in fig. 
3.1 where the asymptotic behaviour has already been subtracted (a similar 
picture appears in |36| ). Because of the zero modes and in accordance with 
Levinson's theorem this subtracted phaseshift starts close to 6tt and falls 
off rapidly on a momentum scale of 2m w . It is essentially this function 



Figure 3.1: Scattering phaseshift from model 
A with asymptotic behaviour subtracted. 



which enters the expression ( |2.26|) for the Casimir energy. The asymptotical 
constants do, at, a<i are known analytically for the N£a model which fact 
served as a test for the program. Numerical values for the full lagrangian are 
close to those given by Moussallam [f36fl . 

Next, we are going to investigate the dependence of the soliton mass in 
1-loop approximation (|2.26| , |2.28|) on the chiral scale /i. In close analogy to 



the 1 + 1 dimensional kink |30] we may expect scale independence only to 
zeroth order in the large expansion parameter of the theory, which in our 
case is N c (an 0{N^ 1 ) scale dependence must in principle be compensated 
by 2-loop contributions). 

In order to accomplish this task we have to choose the scale, say fi = m e , 
such that the lagrangian is dominantly 0{Nc) and then we may look at the 
scale dependence of M(n) in the vicinity of ~ m g . Because we do not know 
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the LECs of ChO > 6 and their scale dependence, we use the pure 4th order 
lagrangian A. 

Before we discuss the numerical results, we will show analytically that the 
mass is indeed scale independent to 0(N®) in the vicinity of \i ~ m e provided 
all ChOs higher than those already contained in the starting lagrangian are 
negligible. 

For this purpose we notice that the LECs of ChO 4 scale according to 
CM = <IW - ^M^) , f i( m s) = 0{N C ) . (3.13) 



In the vicinity of /i ~ m e the soliton mass ( [2.28 ) behaves like 

M„&(aO) ^ M)(^K)) + J-5m4)E7*/^A 

j/7t m g i=1 j 

= M (^K)) + + ^ <n ^) • < 3 - 14 ) 

Because the stability condition for the iVc soliton is calculated from 
Mo(£l(m B )) there is no contribution from the \x dependence of the chiral 
angle. The low energy constants enter the e.o.m for the fluctuations 

and consequently also the phaseshifts 5{p) as well as the asymptotical con- 
stants ratios IM/f 2 which are all of 0{N%). From ( |3~T3| ) it 



is then immediately clear that the \i dependence in all these quantities, E® as 
included, does not appear until 0(N^ 1 ). In ChO 4 the remaining terms in 
( |3.14D just compensate for the scale dependence of the counter terms in ( p. 26 . 
|2.28| ) (a is at least ChO 2) which is then scale independent to 0(N®). 

However the constants m 4 ao, m 2 ai, a<i contain all chiral orders, and the 
higher orders [ChO > 6) are not compensated for. Therefore, the expression 
( |2.26| , |2.28| ) cannot be strictly scale invariant even at (D(Nq). Invariance 
could only be restored if higher ChO terms in the lagrangian were switched 
on upon leaving \i = m g where they were assumed to be zero. This implies 
that the scale dependence calculated numerically does not only comprise of 
small O^N^ 1 ) effects but also measures the magnitude of higher ChO terms 
not accounted for through the usage of an effective Skyrme parameter. All 
the more it comes as a surprise that the soliton mass in tree + 1 loop depicted 
in fig. [O] (a) (solid line) turns out to be almost scale independent over a very 
large region of \i. The scale dependence of the tree contribution (dashed line) 
is nicely compensated for by that of the 1-loop piece. At small scales /i < 
550MeV the onset of scale dependence is rapid till the soliton is destroyed at 
/i ~ 420MeV by the increasing ChO 4 symmetric term (c| = 0.017m^ 2 ,|53|). 
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For comparison, we also evaluated a model with e = 7.24 at [i — 770MeV, 
this choice corresponds to the value implied by the numerical size of the 
LECs at scale /i — m e (table 2.3). 

(a) 



(b) 



Figure 3.2: Scale dependence of the soliton 
mass for models which at scale fi = 770MeV 
correspond to (a) e = 4.25 (g u = 0.0) and (b) 
e = 7.24 (<7 W = 0.0). Dashed lines: tree values, 
solid lines: tree + 1-loop. 
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Any notion of scale independence ceases to exist in this case (fig. |3.2| (b) ) 
which strongly supports the conjecture that the usage of effective parameters 
might emulate missing higher chiral order terms. 



3.2 Scalar properties 

The scalar field which couples to the quark mass matrix is introduced into 
the SU(2) lagrangian (|2.47| ) by the replacement m 2 — > m 2 (l + <;) 

L(e) = L-ef 2 J d\ g[m 2 (l - u) + 2c™u(l - u) - c k A cx^ l u} . (3.15) 

The special choice 

(=io{ ^ tr)= ?^m : (=(p _ pT (3 . 16 ) 



where p and p' denote the nucleon four momenta leads to 

dM(e) 



M(e) = M + ea(t), a(t) 



de 

the scalar formfactor a(t) in the time like region. 



(3.17) 



e=0 



According to the cr-term update [27] the cr-term should lie around a 



cr(0) ~ A5MeV. The interesting quantity for the extrapolation of the cr-term 
from the Cheng-Dashen point is a(2m 2 ) — cr(0) which is estimated in [54 



to be 15.2 ± OAMeV. This value is connected with an extraordinarily large 
scalar radius < r 2 > a = Qa'{t)/a = 1.6/m 2 . 

3.2.1 Tree approximation 

The scalar formfactor in tree approximation is immediately obtained from 
a{t) ^ fl [ d\ z (v / tr)[m 2 (l - c) + 2#c(l - c) + c\{F 12 + 2^)c] . (3.18) 



There are eventually relativistic corrections [|57| which however neither affect 
the cx-term nor the scalar radius 

2 

a fij&r [m 2 (l-c) + 2c™c(l-c) + ^(F' 2 + 2^)c] (3.19) 

^ J d\ r 2 [ml (1 - c) + 2c™c(l - c) + c\{F' 2 + 2 S -)c\ . 



<r 2 tr^e U 
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The main contribution to the a-term is supplied by the familiar pion mass 
term of Ch02 (first term), the Ch04 terms contribute only about lOMeV. 

In the chiral limit m n — > the scalar square radius < r 2 > CT diverges as is 
noticed when the asymptotic behaviour of the chiral angle 



3g A 1 + m n r 

87T/2 r 2 



F ™ - " e~ m - r (3.20) 



is inserted into ( |3.19| ) . In particular, we obtain in this limit 



6a' (0) =< r 2 > a a m *=° [ Sr r 2 F 2 m "=° , (3.21) 

2 J 128tt/^ 



the standard result of HBChPT with external nucleons || |55| . This is due 



to the A states degenerate with the nucleon in soliton models at leading 



order N c 



In general, soliton models, provided realistic parameters are used, overes- 
timate the cr-term and underestimate the scalar radius in tree approximation, 
the latter deficiency being caused by the former. Eventually, loop corrections 
are able to cure both shortcomings. 

3.2.2 Loop corrections to the scalar formfactor 

Here we will discuss the calculation of the Casimir energy in the presence 
of the external scalar field ( |3.16| ). The hedgehog remains solution and the 
stability condition obtained by variation of the tree expression (|3.17|) with 



respect to the chiral angle F picks up contributions from a(t) fl3.18| ) . 
From the terms quadratic in the fluctuations 

L(e) = L-^J d 3 r i (Vir) [-2c k 4 cr) 2 + 2c\c{d i 'n T d i ri T + 8^8^) 



+ 4cl 2C - S 7] L rVr] T + mlcrf - 2c\F' ' s{rf L )' + Ac k A — —7^ 



c k c(F' 2 + + 3r7 2 ) + 2c-(c(2c - l)rj 2 - 2s 2 V 2 ) 



(3.22) 



the e.o.m for magnetic and electric modes may be generated using the stan- 
dard procedure. 



From eq. ( |3.22| ) it is noticed that in the vacuum sector or equivalently 



far away from the soliton centre (F — > 0, s — > 0, c — > 1) the cr-term does 
not vanish. This implies that h^{e) does in general not correspond to a free 
Klein-Gordon equation. However, for t = (zo = 1) and far away from 
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the soliton centre, the sole effect of the terms in ( |3.22| ) stemming from the 
external scalar field is to alter the asymptotical pion mass 



ml( l + e) + 2 c n A 
1 + 2c k d e 



<^K= 'V. =bl 4 (3-23) 



such that hl(e) corresponds again to a free Klein-Gordon equation, but with 
the mass ( |3.23| ) . This change related to the vacuum cr-term 0] is easily taken 
into account by just using the mass (|3.23|) in the expression for the Casimir 
energy ( |2.26| ) . Note here, that the renormalized LECs do not depend on the 
external field which simplifies the situation considerably. 

The external scalar field ( |3.16| ) is rotationally invariant and, for t — 0, 
translationally invariant, too. Therefore all zero modes must be recovered in 
the presence of the external scalar field with t — at zero energy. This serves 
as a crucial test that we have solved the stability condition for F and the 
e.o.m for the fluctuations correctly. 

For i ^ 0, i.e. for the scalar radius (io(\/tr) ~ 1 + r 2 t/6) the terms 
proportional to the external scalar field (|3.16|) lead to a confining potential 
(oscillator potential in the case of the scalar radius) and we are not able to 
use the phaseshift formula ( (2.21|) without manipulation. In coordinate space 
the expression / Sr (h(e) — ho(e)) collects contributions only from a region 
r < R where the soliton profile is unequal zero. We integrate the phase 
shifts 5(e) and S (e) of h 2 (e) and hl(e), respectively, to a radius R where the 
asymptotical mass becomes 

m 2 = ml(l+et (ViR)) + 2c?et (ViR) 

m ° l + 2c k 4 ei {ViR) ' 1 ' ' 

The difference 5(e) — 5 (e) entering the phase shift formula (|2.21|) should 
then lead to a Casimir energy independent of R, provided R was chosen 
large enough. 



3.3 Axial properties 

The most general form of the nucleon axial current in the Breit frame p p 



(E, —\q), V w — {E, \q), E = y M 2 + \q 2 (no energy transfer = (0,q) ) 
is given by 

(N r (lq)\A a \N(-±q)) = {N'^G^q^q^N) 
(N'(\q)\A<t\N(- l -q)) = (N'\ ( ^G A (q 2 ) (a, - q % (qcr)) (3.25) 



HG A (q 2 ) ~ -^ r M P (q 2 ))q l (q ( r) ±\N) . 



2 



AM 2 



Tn 



2 
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The time component Aq and the corresponding pseudotensor formfactor 
Gt{q 2 ) are related to that part of the axial current which has positive G- 
parity (so-called second class current). In soliton models the time component 
of the axial current in tree approximation 

A a tree + ^ sc ( f X 

f 2 1 

= ^(br + c{o) sc Spigh-iDap, R q } (3.26) 

is proportional to the angular velocity Q R and to the right angular momen- 
tum R, respectively. Because after hermitean ordering the operator 

(N'\e pig ^{D ap ,R g }\N) = l -{N'\R 2 D ai -D ai R 2 \N) = (3.27) 

vanishes between nucleon states there is no contribution to the axial pseudo 
tensor formfactor ||18|| . There may well be nonvanishing loop corrections; 
however, loop corrections to nonadiabatic quantities (as (|3.26|) ) become very 
much involved and are not treated in this paper. 

The O(Nc) spatial components of the axial current appear in soliton mod- 
els as A\ = D av A\ where the Euler angle dependent D-function transforms 
the intrinsic A\ to the lab. system. For the nucleon matrix element we may 
use D ap = —\T a a p and the decomposition 

A\ = A L hr p + A T {5 ip - fif p ) (3.28) 
into longitudinal and transversal parts to obtain 

(N'\ J d\ e^A^N) = (N'\ J d\ e iC L r D a pA p AN) 

= ~{N'\ j d\ e^lA^rp 

+A T (5 ip - rirp)]cr p -^-\N) 

= ~(N'\ J d\ [ Jo (qr)(A L + 2A T )a l 

+ j2 (qr)(A L - A T )(a t - 3ft($<r))] j\N) (3.29) 

The decomposition ( |3.28|) is evident for the tree approximation (see ( [3.33j ) 
below) and holds still for the 1-loop part when adiabatic fluctuations depend- 
ing on r and momentum p are inserted and the momentum is integrated out. 
Then A\ is again a function of only r and the decomposition ( |3.28D is possible. 

Comparing (|3.25|) and ( |3.29|) the axial formfactor Ga and the induced 
formfactor Gp can be read off as 

E 2 { 

— G A {q 2 ) = -gJ d 3 r [(jo(qr) + j 2 (qr))A L + {2 Jo {qr) - j 2 {qr))A T ] 
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,2 



G A (q 2 ) ~ ^f 2 Gp(q 2 ) = -g / d\ [(jo(qr) - 2j 2 (qr))A L (3.30) 

+2(j (qr) + j 2 (qr))A T ] . 

Because Gp is not very well known experimentally we will concentrate on 
Ga, in particular on the axial vector coupling constant 

9a = G A (0) = ~j d\ (A L + 2A T ) = ~\j d\ A\ . (3.31) 

Relativistic corrections |j7] may be accounted for through the replacement 

GM) - ^Gx(^). (3.32) 

For small momentum transfers these corrections are of minor importance; qa 
is not affected at all and the axial square radius experiences only a slight 
increase by 3/4M 2 = 0.03/m 2 (compare (|3.36|) ). 



3.3.1 Tree approximation 

In tree approximation, the spatial components of the intrinsic axial current 



£P tree p 



SC 



(3.33) 



(b L + c{u)F'fir P + (6r + c{U) — (5 ip - f % f v 

r 

are of the form (|3.28|) as already mentioned and the axial formfactor is ob- 
tained as 

E. Ga ^) ^JJl J Sr [tio + j2)(b L + c{q 2 )F' + (2j -j 2 )(b T + c{q 2 )^} . 

(3^34) 

For completeness, we give the corresponding formulae for the axial vector 
coupling constant and the axial square radius 



9a 



tree 

9 



^ Jl \ *™ e ^ ^ f 7T 



( d\ (b L F' + 2b T —) (3.35) 
J r 

J d 3 rr 2 (b L F' + 4b T ^) 



4M 2 !5g A 



+c{A U\ (F' + 2-) . (3.36) 
' J r 



%9a 

It should be mentioned that the main contribution to the axial radius is due 
to the non- minimal coupling (last term in ( [3.361 ) ) . In tree approximation the 
axial coupling constant turns out too small in chiral soliton models provided 
reasonable parameters are used and this deficiency can not be cured through 
inclusion of vector mesons. Therefore this quantity is considered to be one 
important candidate for which loop corrections may prove to be essential. 
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3.3.2 Loop corrections to the axial formfactor 



Here we want to consider the loop corrections to the axial formfactor ( |3.34j ). 
For this purpose we have to discuss the Casimir energy in the presence of an 
external axial field af 



L{e) = L + e Jd 3 r a^A^ = L + e J d\ afl?, (3.37) 



A\ — DapA?, a" — D ap a? 



which finally leads to the axial formfactor 5u(<? 2 ) = E/MGa{q 2 ) 

dM{e) 



M(e) = M + eg A , g A (q 2 



ds 



(3.38) 



e=0 

According to ( |3.28| ) and ( |3.30| ) the external axial field has to be chosen as 

2 1 

5 * = g [O'o + h)nr p + (jo - 2^2) [S ip - r l f p )] 

2 1 

a i = g [(Jo + j2)Dapfif p + (jo - ^32){Dai - A^tfj,)] • (3.39) 

This simplifies for the axial vector coupling constant g A = S'a(O) to 

2? = ^p, < = \Dai. (3.40) 

We wrote down here the axial currents and fields in the lab frame (without 
tilde) just to display their correct axial transformation properties although 
we do not need them in the following: the whole calculation is performed in 
the intrinsic system using the standard procedure. 

If we expand the lagrangian (|3.37|) in the presence of the external axial field 
to linear order in the fluctations, we notice that again the hedgehog ansatz 
solves the static e.o.m . The stability condition is obtained by variation of 
with the axial formfactor in tree approximation (|3.34p which yields a 



g-dependent chiral angle and makes the linear term in the fluctuations vanish. 

In close correspondence with the longitudinal and transversal operators 
appearing in the external axial field we can distinguish two pieces in the 
lagrangian quadratic in the fluctuations 



L{e) t L - 2 ^{jSr{j Q +j 2 )[F>rf T + - 



—r) L rVr) T + — (2r/£ + rfr) + 



(3.4i; 
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from which the e.o.m for the fluctuations may be generated (for the sake 
of simplicity we have listed only terms produced by the N£cx model). From 
eq. (|3.41| ) it is clear that the axial current matrix element is zero in the vacuum 
sector F — > 0. Therefore h^(e) = corresponds to a free Klein-Gordon 
equation with pion mass m n and we do not have the difficulties discussed in 
the preceding section for the cr-term. 

For the 1 loop calculation of g& there exists an alternative method con- 



cerning the evaluation of the integral (|3.37| ) which we shall describe briefly. 
The divergence of the axial current is related to the symmetry breaking terms 
in the lagrangian fl2.47j ) 



0M£ 



An + 0M? 



mlD. 



ap 



fl s ^ P + U(cr p r] L + 7] pT ) - -f p rf 



+ 



(3.42) 



where we expanded the standard pion mass term up to quadratic order in 
the fluctuations (similarly for the other symmetry breakers not listed here). 
We may now use eq.( |3.42| ) at order rf to simplify the integral in ( ^.37| ) by 
partial integration, 

L(e) t L + ^ / d\ 8%!? 
9 J 

= L-j J d\ Xpd'A? 

= L + ^ f d\ rsrf + ■■■ , (3.43) 
9 J 

where we have checked that the surface terms vanish. Because we neglect the 
contribution from Aq which is 0(Nq 1 ) the loop corrections to calculated 
according to ( |3.41| ) and ( |3.43j ) , respectively, differ just by ^(A^^ 1 ). 

The disturbing result is now that this difference turns out to be large in 
our soliton models (.45 for parameter set A) which suggests that the O^N^ 1 ) 
contribution to gA might be large, too. The consequence of a large 0{N^ 1 ) 
contribution to gA and it's connection to current algebra will be discussed in 
the following subsection. 



3.3.3 Current algebra and l/iV^-expansion 

From the time component of the axial current in the infinite momentum 
frame and the chiral charge commutator 

[QlQ b 5 ] = ie abc L c , (3.44) 



45 



Kirchbach and Riska derived a model independent version of the Adler 
- Weisberger sum rule 

g\ = l + R. (3.45) 



In ||58fl , the "1", produced by the commutation relations, is of 0(N C ) and R 
is always positive. 

From ( |3.45| ), it is immediately obvious that the 1/Nc expansion cannot 
converge reasonably for both, g 2 A and R, because otherwise we would expect 
the 0(N C ) contributions to both of these quantities to be small compared 
to the experimental values of 1.56 and .56 respectively, which violates ( |3.45| ) 
at 0(N C )\ Obviously, the 1/Nc expansion can converge rapidly, if it is to 
converge at all, only for one of the two quantities. From our soliton model 
calculations we conjecture that this is the case for R and not for g 2 A : For 
example, the model A in tree + 1-loop provides the 0(N 2 ,) and O(Nq) 
contributions g\ = R= .83 — .41 = .42. Clearly, this seems to converge to the 
experimental value of R rather than to that of g\. Interestingly, such a result 
(negative 0(N C ) contribution to gA) was already found in |TL|, although we 
disagree with the procedure used there. If we add the current algebra "1" 
of 0(N C ) from eq.( |3.45| ) by hand, thereby assuming that gA calculated up 
to 0{1/Nc) exhausts relation (|3.45|) , we obtain the quite satisfactory result 
gA = -91 — .25 + .54 = 1.20 at the expense of a large 0(N^ 1 ) contribution to 
gA- The suggestion that there might be a sizeable 0{1/N C ) piece contained 
in gA is supported by the alternative calculation in the preceeding subsection, 
which yields a large difference precisely in this order. However, an explicit 
calculation of this piece, which involves tree (with two angular velocities), 
nonadiabatic 1-loop (one angular velocity) and adiabatic 2-loop contributions 
seems forbiddingly complicated. Using instead the ad hoc addition of the CA 
" 1", we may estimate the contributions to for the other parameter sets, 
too. The individual 1/Nc contributions to gA for the various models are 
listed in the table |3TT]. It is noticed that the 1/Nc piece is large and positive 
throughout and increases the tree + 1 loop value towards the experimental 
datum, except for the case of parameter combination D (e = 3.75, g u = 
0), designed to yield acceptable numbers for already at tree + 1 loop 
level at the expense of all other quantities, which turns out to considerably 
overestimate gA and should therefore be discarded. 

For the axial radius we may employ the same procedure. Instead of ( 3.45|) , 
one has p9 | 

g\ < r 2 > A =< r 2 >\ +X (3.46) 

with the isovector square radius < r 2 >\= -58fm 2 which determines the 
nucleon-A transition radius < r 2 >* A . (Note that, in the end, we will be 
interested in the quantity gA < t 2 >a which is the slope of the formfactor). 
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The isovector square radius related to the isovector charge square radius and 
magnetic moment 

<r 2 >r=<r 2 >^-A ( ^_l) (3.47) 

appears due to the CA commutation relations and is 0(N C ). 

From the experimental datum < r 2 >a= Alfm 2 we again find a small 
value [gA < r 2 >aY 1+0 ^ = -19/m 2 (superscripts here denote N c orders) in 
tree + 1 loop and a large positive 1/N C correction [g A < r 2 >a] < ~ 1 - ) = -34/m 2 . 
Again for the reasonable sets (A) — (C) the pattern is repeated and the tree + 
1 loop result enhanced toward the experimental value gA < t 2 >a= -53/m 2 
although in this case there is an overestimation for the first three models 
and model (D) is not widely off the mark. Altogether, for the axial radius, 
numbers are less conclusive in view of the large error of the axial radius of 



Table 3.1: Comparative listing of tree, 1-loop and 
estimated 1/Nq piece of g& and < r 2 >a for the 
four parameter combinations A, B, C and D. 





A 


B 


C 


D 


(1) 

9a 


.91 


.96 


1.0 


1.14 




-.25 


-.15 


-.11 


.14 


si" 11 


.54 


.48 


.44 


.34 


9a 


1.20 


1.29 


1.36 


1.62 


[9A<r 2 > A r } [fm 2 ] 


.41 


.42 


.41 


.62 


[gA<r 2 >A] [0) [fm 2 ] 




-.13 


-.07 


.11 


-.16 


[g A <r 2 > A ] (_1) [fm 2 




.38 


.35 


.28 


.28 


[g A <r 2 > A ] [fm 2 ] 


.62 


.65 


.76 


.69 


<r 2 > A [fm 2 ] 


.55 


.54 


.59 


.46 



This unusual scenario (negative O(Nq) and large positive 0{N^ 1 ) contri- 
bution) for axial quantities is in striking contrast to all other models (nonrel- 
ativistic quark model, bag models, NJL model) which find a positive 0(Nq) 
and no O^N^ 1 ) correction. This is very serious because not only the exper- 
imental value of gA itself but also the individual contributions of different 
Nc orders are fixed in principle and if the soliton approach is indeed cor- 
rect, all others are wrong and vice versa. Interestingly, the same numbers for 
the quantum correction at 0{N^ 1 ) come out of a calculation which assumes 



47 



that Skyrme type models should give a positive remainder R and imposes 
SU(4) current algebra for the spatial components of the axial current. This 
calculation is presented in some detail in appendix A. 

As an addendum, we have calculated the scale dependence of qa starting 
from set A. The result, depicted in fig. [O] shows the scale independence 
to be not worse than for the mass considering the fact that the tree value 
(dashed line) of g& is much more sensitive to the change of scale than the 
soliton mass. Due to the too large symmetric 4th order term which incurs 
numerical difficulties, we were unable to calculate the correction at scales 
below QOOMeV. 



Figure 3.3: Scale dependence of g& for a model 
which at scale \i = 770MeV corresponds to 
e = 4.25 (g u = 0.0). Dashed line: Tree value, 
solid line: tree + 1-loop. 



3.4 Electromagnetic formfactors 



The nucleon matrix element of the electromagnetic current in the Breit frame 

(3.48) 



(N'^q)\M0)\N(~q)) 



(N'(lq)\Ji(0)\N(~q)) 



(N'\G v E (q 2 )rs + G%(q')\N) 
G v M (q 2 ) , GlM 2 )' 

T 3 



(N'\ 



2M 



2M 



i((T X q)i\N) 
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fixes the isoscalar and isovector formfactors which are linear combinations of 
the proton and neutron formfactors 



qV^B.M "t" ^E,M) ) ^ E,M 



2\^E,M ^E,M) 



(3.49) 



normalized to the isoscalar and isovector charges and magnetic moments 

(3.50) 

The corresponding radii 

6 dGZM 2 ) 



< r 2 > S ' V - 



O sy (a 2 



dq 2 



(3.51) 



q 2 =0 



are related to the slopes of these formfactors. 

According to the charge operator Q = ~(r 3 + I /No) which guarantees 
integer charges of O(N^) for mesons and baryons the electromagnetic current 
in soliton models 

J, = V* + \v° = D 3p V? + iy M (3.52) 

decomposes into an isovector part represented by the third component of the 
vector current and an isoscalar part related to the baryon current. For the 
nucleon matrix element we may again use D 3p = — ^r 3 a p to obtain 

(N'\ J d\ e^ r J,\N) = (N'\ J d\ e^ r [~T 3 a p V* + \v° 

Comparison with ( |3.48| ) gives the electromagnetic formfactors 

G V E (q 2 ) = ~\j d 3 r j (qr)a p V p 



\N) . (3.53) 



GW) 

GW) 
2M 

G S M (Q 2 ) 
2M 



d\ Uqr)V ° 



- 1 d \ 



qr 
jijqr) 
qr 



x 



(a x r)^ c 



(3.54) 



Relativistic corrections [[57]] may be taken into account by the replacements 



G S /{^ 



G 



s,v 



M 2 q 2 



, G S m (q 2 ^ 



M 2 ^ sy ,M 2 q 2 



r 1 ' 



(3.55) 



Again, for small momentum transfers these corrections are of minor impor- 
tance; charges and magnetic moments as well as the electric radii are not 
affected, the magnetic square radii collect a small contribution 3/2M 2 = 
0.06/m 2 . 
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3.4.1 Tree approximation 

In tree approximation we obtain for the intrinsic isovector current 
VS tr = ^(b T + c{n)s 2 {R p -f p (fR)} 



V[ = e fAb T + dnf-e^rj (3.56) 



r 

and similarly for the isoscalar current 

1 1 9F'e 2 

mf, Aix A r z 

UJ 

V? '= (l-\n)B it B^^rxR)^ (3.57) 

where the nonminimal ChO 6 coupling through the a;- meson, essential for 
the isoscalar radius, was taken into account. With the currents (|3.56| , p. 57] ) 
the formfactors ( |3.54 ) read in tree approximation (R p = —a p /2 for nucleon 
states) 



2 



G V E (q 2 ) = ^1^ Jo(qr)(b T + c{q 2 )s 2 G V E (0) 

G S E {q 2 ) = \{l ~^)J d\ j (qr)B Gf (0) = \ (3.58) 

UJ 

G v M {q 2 ) = 2M& J d\ 3 -^(b T + c{q 2 )s 2 GlM = f = f 9 

GW) = £(1 - 4) ^r 2 B <&(0) = ^ = ^<r 2 > B . 



66 v m*' J qr u JWV ' ^ 60 

The corresponding radii are 

6 



< r z >^ = < r z > e _c 4-q- Jars, <r > e = — y a r r 6 T s 

< r n > B = J d 3 r r n B 



< r 2 >| = < r 2 > B +— , 



< r" > 



29 3 3 < r 4 > B 6 .„ „. 

<^ >m = 777^ + ^— T— + (3-59) 



2M 2 5 < r 2 > B mj 

where the relativistic correction due to ( |3.55| ) was included. The nonminimal 
couplings through the ^-meson in ChO 4 and through the u- meson in ChO 6 
are essential and lead to a considerable enhancement of all these radii. 
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In the following we will calculate loop corrections to some of these quanti- 
ties. As is noticed from Q3.56] , |3.57|) the adiabatic ones containing no angular 
velocity repectively no angular momentum for which the method to calcu- 
late the Casimir energy outlined above applies are the isovector magnetic 
formfactor of 0(Nq) and the isoscalar electric formfactor of O(Nq). 



3.4.2 Loop corrections to the isovector magnetic form- 
factor 

For the calculation of the Casimir energy in the presence of an external 
isovector field vf 

L{e) = L + eJd 3 r v?V? = L + e J d\ dfVf, (3.60) 

V a = D V p v a = D 77? 

we proceed in very much the same way as for the axial field of the previous 
section. The external field which leads to the magnetic formfactor <?J^((/ 2 ) = 
G v M (q 2 )/2M in the form 



v v . 2 , dM{e) 



M(e) = M + eg v M} g v M (q 2 ) <y 



(3.61) 



e=0 



follows from flpp and ( $M ): 



v 



7: £iptxe, v, t = - D 3p s ip ixi. (3.62) 

2 qr 2 qr 



The isovector field vf in the lab. frame has the correct transformation prop- 
erty. Again the hedgehog solves the static e.o.m. and makes the term linear 
in the fluctuations vanish. From the term quadratic in the fluctuations, 



L{e) t L + ^Jd 3 r^[(c 2 ~s 2 )(2 V 2 L + r 1 2 T ) + 2rc VL Vr lT + 



£ f ^ 3ii 
qr 

(3.63) 



the e.o.m. for the fluctuations are set up. Here, apparently the additional 
terms due to the interaction do not vanish in the vacuum sector such that 
hl(e) = h"l + ew with an asymptotic potential w ^ 0. However, for the 
Casimir energy we have to evaluate / dfr (h(e) — ho{e)) — J d^r (h(e) — ho — 
ewo/2h ) = J d\ (h(e) — h ). The term which contains the asymptotic poten- 
tial turns out to vanish because of / d\ (swo/2ho) = 0, which implies that 
the vacuum matrix element of the isovector current is zero, as it should be. 
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For practical reasons, we must subtract the corresponding term which is also 
zero in the soliton sector from h 2 (e) in order to obtain stable phase shifts. 

In the presence of the vector field fl3.62| ) which is neither rotationally 
nor translationally invariant (even at q = 0) the e.o.m. for the fluctuations 
unfortunately cannot be checked by the zero modes. 



3.4.3 Loop corrections to the isoscalar formfactor 

The external isoscalar field v® = ejo(qr), 

L(e) = L- £ -jd\ vX° = L - |(1 - i) / d\ v° B 

UJ 

immediately leads to the isoscalar electric formfactor 



(3.64) 



M(e) = M + eG%, G%(q 2 



3.65) 



£=0 



Of course, the hedgehog is again solution and the e.o.m. for the fluctuations 
are obtained from 



(3.66) 



where we have used partial integration. It is noticed that for q = the loop 
contributions in ( |3.66| ) vanish identically. This is of course to be expected 
because the isoscalar charge G S E (0) = 1/2 should not be subject to quantum 
corrections. Apart from this (trivial) test we checked the e.o.m. by the rota- 
tional zero mode, which remains at zero energy because the external field is 
rotationally invariant. 



3.5 Polarizabilities 

The electromagnetic polarizabilities of the nucleon have been calculated in 



various approaches, amongst them nonrelativistic quark- |p0|-||62|, MIT bag- 
JkJ and chiral quark models ^§ as well as BChPT § @ [§§ 0. The 
static polarizability is a measure for how easily an electric (magnetic) dipole 
moment may be induced by a constant external electric field e (magnetic 
field b). In the lagrangian this leads to additional terms 

L{e, b) = L + \(A-na h )e 2 + \{A-K(3 b )b 2 (3.67) 
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Experimentally there is information about the dynamic polarizabilities ob- 
tained from Compton scattering |69| , which have to be corrected for 



relativistic and retardation effects caused by the finite size of the nucleon |70 
73] and about the static electric polarizability of the neutron from neutron 
nucleus scattering JF3[ . There is no direct measurement of the static polariz- 
abilities of the proton. The experimental results for the neutron and proton 
electric polarizabilities a n = 12.0±3.5- 10~ 4 /m 3 and a p = 7.0±3.5-10~ 4 /m 3 
are arranged into the linear combinations 

a = ^{a n + a p ) = 9.5 ± 5.0 ■ 1(T 4 /W 5 

m = a n - a p = 5.0 ± 5.0 • lO -4 /™ 3 (3.68) 

which turn out to have different Nq behaviour in chiral soliton models. In 
these models the adiabatic quantity a is leading O(N^) in tree approxima- 
tion (section 3.5.1) and quantum corrections of 0{N^) (section 3.5.2) may 
be important. The non-adiabatic quantity a« comes with an angular velocity 
in tree approximation [73J and gives rise to a 0(Nq 1 ) splitting of neutron 
and proton electric polarizabilities (section 3.5.3). For the magnetic polariz- 
abilities this is quite analogous, only the experimental information is much 
poorer there. All previous attempts to attack the problem in the framework of 



chiral soliton models |20| [^TJ |E| HU |25j rer y on the adiabatic tree approx- 
imation. Trivially there is no splitting of the neutron's and proton's electric 
polarizabilities (ao = a/3 = 0). 

More severely the numerical values obtained, in purely pseudoscalar mod- 
els a ~ 18 — 30 • 10 _4 /m 3 depending on the parameters used, come out much 
too large. Nonminimal couplings to the photon do not help here, they are of 
negligible influence (section 3.5.1). This finding is at variance with the seem- 
ingly much smaller results a = 10.4 — 14.8 • 10 _4 /m 3 obtained by adding 
vector mesons explicitly which include nonminimal couplings naturally. 
The reason for this discrepancy is due to Qq an d components in the vector 
meson soliton profiles induced by the electric field e which lead to additional 
O(Nc) contributions to the electric polarizability and which are simply ne- 
glected in the approach mentioned |2"2]| . The situation is easily understood 
and the missing terms can be traced by looking at the local approximation 
which converts the vector meson into a purly pseudoscalar model. Taking 
these contributions properly into account, vector meson models would end 
up with large electric polarizabilities comparable to those of the correspond- 
ing pseudoscalar models. 

For the magnetic polarizabilities the situation is even worse. Numerical 
values obtained so far range from (3 = —7 — (+8.5) • 10 _4 /m 3 |^TJ [23[] [24 



Up to now there does not even exist a correct and complete calculation of 
the leading O(Nc) contribution in the simplest Skyrme model. 
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3.5.1 Electric polarizability in tree approximation 



The isovector photon field Vq = —er leads to a contribution in the lagrangian 
quadratic in the electric field e which in tree approximation is 0(Nq) 

l'2 2 

L~ r/ = J -^- j d\ {er) 2 s 2 b T D 3a D 3b {5 ab -f a f b ) 

+ fl{2c\ - c{)eV J d\ s 2 D 3a D 3b (5 ab - r a r b ) (3.69) 



(We use f^r = t|? throughout). There is no such contribution for the isoscalar 
photon field and most importantly there is no O(N^) term at all which 
is linear in the electric field as e.g. in vector meson models which would 
induce adiabatic soliton components driven by e. There is, however, such a 
term in O(Nq), nonadiabatic, which contains an angular velocity and which 
consequently splits neutron and proton polarizabilities. This term is treated 
in subsection 3.4.3. Concluding, in leading 0(Nq) eq. ( |3.69|) represents the 
complete contribution to the electric polarizability 

f2 2 . f2 2 



187T 



/ d\ r 2 s 2 b T + ^-(2cl - c{) / d\ s 2 . (3.70) 
J on J 



Here we used that for nucleon states D 3a D 3b = ^5 ab . Alternatively we could 
have averaged over the orientations of the electric field e which leads to 
the same result. For the calculation of the quantum corrections in the next 
subsection we will have to do both simultaneously, taking nucleon matrix 
elements and averaging over the e-orientations. 

The second term in (|3.70|) represents the non-minimal couplings. If the 



chiral scale is chosen such that the LECs t & ~ 2£ r 5 ~ —l/8g 2 are exhausted 
by the ^-meson resonance E7| then the contribution from the non-minimal 



couplings c{ ~ 2c% vanishes exactly. The Gasser Leutwyler LECs at \i — m 
lead to an unsignificant enhancement of the electric polarizability due to the 
non-minimal couplings. Thus, in contrast to the common belief, neither non- 
minimal couplings nor the explicit inclusion of vector mesons is able to cure 
the problem with the much too large electric polarizabilities obtained in tree 
approximation from soliton models. 

Such a conclusion is of course valid only if the way to calculate polariz- 
abilities as sketched above is correct, which has recently been doubted ]75 



76]. The explicit statement was that the terms quadratic in the external elec- 
tric or magnetic fields (seagull terms) do not contribute at all. Since seagull 
terms are in fact the only contribution at leading order Nc, this claim, if 
true, would make quite a difference, and we shall in the following discuss it 
in some detail. 
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A first hint that not all is well with the statement voiced in [[HJ [|7fJ may 
be derived from the fact that the tree expression for the electric polarizability 
( |3.70D evaluated in the chiral limit (where an analytical treatment akin to the 
case of the a- term is possible) agrees with the result obtained in HBChPT 



a 



9 



dr r A F 2 



5g 



(3-71) 



4vr 32 7 r/> 7r ' 

In order to put our concerns on a more solid footing, we insert the ansatz 
(|2.14 ) for the time dependent matrix U into the lagrangian coupled to the 
external field Vq (intrinsic system). It's relevant terms then read 



L(e) = -M + - J d\ (v a n 2 ab f) b - Vah 2 ab r] b ) 



+ 
+ 



ie/dV(nf + ev5Kn^(n? + 



ev. 



(3.72) 



with the moment of inertia and z e a the normalized rotational zero mode 

(PD, 

/l - (3.73) 



^-a.p.iTi 



The last term explicitly written down in p. 72] ) contains the rotational energy 
as well as the seagull contribution. The collective rotation (translations are 
unimportant in the present context) together with the fluctuations rj gives 
rise to a redundancy which shows up in the conjugate momenta 



Tin 



5L_ 
Srja 



n 



ab 



Vb 



(3.74) 



Re 



dL 



e / d\ zy ab rj b - e / rfv zy ab z{(nf + e%) 



making them linearly dependent (so called primary constraint [^0| 

R e = V® j ' d\ Z e a 7T a . 



Due to this fact, the naive Legendre transformation (used in |75 
hamiltonian 



(3.75) 
to the 



H = M + 



1 



d s r 



7r a n ab 7T b + r] a h 2 ab r] b ] + eV& / d\ v^z[ir a (3.76) 



is not well behaved: The collective coordinates have been removed completely, 
and gone are the seagull term as well as the soliton's rotational energy. 
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In order to keep the collective coordinates, one has to impose secondary 
constraints |40| conjugate to the primary ones ([3.75 ). There is an obvious 
choice for these: the fluctuations should be orthogonal to the zero modes, 
which have already been accounted for through the collective variables 



d 3 r zC a n lbVb 
d 3 r 4 n lbVb 







d\ z c „-K n = 0, 



Tin 



n lbVb 



(3.77) 



These constraints have to be added with multipliers to the lagrangian and, 
going through Dirac's procedure f40| , one finally obtains the proper hamil- 
tonian. However, the same objective may be achieved faster by decomposing 
the conjugate momenta (|3.74|) 



Tin 



7T„ 



Tin 



717, 



Qnl b z[(nf + ev f ) 



(3.78) 



into collective and fluctuational parts and inserting them into the hamiltonian 
( |3.76| ) . It is noticed that the collective motion decouples completely from the 
fluctuations 



H 

H 



H. 



n 



ii 



M + 



20 



<rr i^an^ 7T b + r] a h ab 7] b ) 



E 1 = eR e J d\ z e a n 2 ab z{v{ 



2 



d3r z a n lb4vo 



(3.79) 



d 3r zAdvlvl 



The rotational energy R 2 /2Q is recovered as well as the seagull term i7 77 . So 
far we have not specialized on the particular photon field Vq. For a spatially 
constant potential Vq = const., if 77 does indeed vanish, as is to be expected 
since Vq = const, corresponds to a global rotation. For the constant electric 
field Vq = —erDza needed to calculate the electric polarizability the first 
integral in if 77 vanishes for parity reasons and we end up with H 1 = and 



H, 



J Sr {er) 2 s 2 b T D 3a D 3b {5 ab - f a f b ) , (3.80) 



the familiar soliton model result, compare eq.( ^.69[ ). Therefore we definitely 
do not agree with the statement that the seagull terms should vanish and 
are confident that the result (|3.70|) for the tree contribution to the electric 
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polarizabiliy in soliton models is indeed correct. Similar conclusions were 
independently drawn by Scoccola and Cohen |77j 

In the next section we will provide the method to calculate the quantum 
corrections to this interesting quantity in order to learn whether these 1-loop 
contributions are eventually able to reduce the too large tree level values. 



3.5.2 Loop corrections to the electric polarizability 

The calculation of loop corrections to the electric polarizability is very much 
involved, concerning both, the theoretical input and the computational effort. 
Although the output is only one number we decided to undergo this pain, 
because this quantity persistently comes out too large in tree regardless of 
the model considered and therefore may serve as a crucial test of the loop 
expansion. 

In the end we want to calculate our baryon property from the Casimir 
energy 

1 dM(e 2 ) 



M(e 2 ) = M - 2ne 2 a, 



a 



2tt de 2 



(3.81; 



in the familiar way. For this purpose the photon field v\ = —er has to be 
coupled to the lagrangian (via the covariant derivatives) which then ought to 
be expanded up to quadratic order in the adiabatic fluctuations. Of course 
if we consider a certain fixed direction of the electric field the soliton gets 
deformed under its influence and our whole concept, which is based on the 
spherical hedgehog solution would fail. Luckily we are interested only in the 
static electric polarizability which, defined as the coefficient of e 2 (|3.67|) , is 



independent of the direction of e and we are allowed to average over these 
directions. In addition we may use the relation D 3a D 3b = ^5 ab for nucleon 
states. With these two allowed manipulations we are able to show that the 
hedgehog still solves the e.o.m. in the presence of the electric field e and makes 
the term linear in the fluctuation vanish if the stability condition for the chiral 
angle is calculated by variation of the tree contribution ( |3.70| ) contained in 
( |3.81| ). So far everything is quite similar to the quantities calculated in the 



previous sections. The complications arise with the terms quadratic in the 
fluctuations, or more precisely, with the e.o.m. for the fluctuations. 

The lagrangian expanded to second order in the fluctuations 
L 7 = eD 3a Jd 3 r (er)[(rj X rj) a - 2(1 - c)-q L {r X r)) a + ---\ 

2 \ f r 

£ 77 - -e 2 D 3a D 3b J Sr (er) 2 \r) 2 5 ab - r) a r] b - s 2 {2i] 2 L + rf T )(5 ab - f a f b ) 
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-(1 - c)rj L {f a rj hT + r] aT r b ) + 



(3.82) 



contains a piece linear in the electric field with one time derivative and a piece 
quadratic in the electric field (for simplicity again only the N£a contributions 
are listed). The second term can be simplified by using D 3a D 3b = ^5 ab for 
nucleon states and averaging over the electric field orientations (sr) 2 — > 

l F 2 2 

3 

Naively one might think that the first term vanishes by averaging over the 
electric field orientations but this is not the case for the e.o.m. as we will 
demonstrate in the following. 

The lagrangian ( |3.82 ) introduces additional terms ~ e and ~ e 2 into the 
e.o.m. which may be written formally 



h 2 i] + 2iew 1 fj — e 2 w 11 rj = —n 2 {e)fj 
with hermitean differential operators u> 7 and u> 77 
ieD 3p £r [e pab 

hp^3q 



(3.83) 



(lf 77 )aft 



e 2 D 3n D 3n (£r) 2 



(1 - c)f q (f a e pqb - E pq ah) H ] (3.84) 

S pq S ab - 5 ap 5 bq - s 2 (5 pq - r p f q )(5 ab + f a f b ) 
+ (1 - c)(r a r p 5 qb + 5 ap f q r b - 2r a r p r q r b ) H ] . 



Differentiations appear for the higher chiral order terms omitted in ( |3.82j ). For 
w 7 alone we obtain 9 terms of different isospin structure from the lagrangian 
( |3.67D .The single time derivative in eq.( |3.83j ) is exactly of the form discussed 
in eqs.( p.33| , |2.34|) and consequently 



d 4 x{£nQ) 



T 
2" 



(fr (Jh 2 - e 2 (w 77 - w 2 )) 



(3.85) 



there is no term linear in the electric field e. Thus considering the trace log 
( |3.85| ) we may equally well solve the e.o.m. 



h 2 r] — e 2 {w 11 — w^n 2 w 1 )rj 



-n 2 T) 



(3.86) 



instead of ( |3.83| ). This is of the standard form with an additional potential, 
simple only in the case of the N£cr model (n = 1) 



e 2 D 3p D 3q (er) 2 



Spqdab — 5 ap 5 b 

s 2 S pq r a f b + (1 - c)(r a r p 5 qb + 5 ap f q f b 



IpqVab VapVbq & ^i^j^-pia^-qjb 
-2) 



-(1 - cff a r p r q r b ) + • 



(3.87) 



In general we had to work out, according to ( |3.85| ) and the longitudinal and 
transversal metric, 2x9x9 terms. From Q3.85| ) and (|3.87|) it is recognised 
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that in contrast to the individual potentials U7 77 and w^n 2 w 1 their difference 
is well behaved for r — > oo (F — > 0). 

Finally we may simplify the potentials using again D 3p D 3q = \$ pq and 
averaging over e-field orientations to obtain 



[w 



77 )ab 



w 



7 jab 



— r 
9 



— r 
9 



25 ab - 2s 2 (5 ab + f a f b ) + 
2S ab - s 2 (8 a b + r a r b ) + ■ 



(3.88) 



With w 77 alone we would end up with an infinite loop correction to the 
polarizability of the vacuum which of course should be zero. This is how 
we noticed that we need the first term in (|3.82|) linear in the electric field 
in order to get a reasonable finite result for the quantum correction to the 
electric polarizability in the vacuum as well as in the soliton sector. 



3.5.3 Neutron-proton split of the electric polarizability 

In this subsection we present a derivation of the neutron-proton split a« of 
the electric polarizability in chiral soliton models. A similar approach for a 
chiral soliton model with quarks was published by Broniowski et al. [f74|] . 

The quantity under investigation must contain an angular velocity in order 
to split the values for neutron and proton and is therefore of lower order Nc 
compared to the tree and 1-loop contributions to the electric polarizability 
discussed in the previous sections. 

We consider the term linear in the electric field (non-minimal couplings 
do not contribute here) 

L 7 = ej d\ (er)(D 3a V a + ^B ) (3.89) 

which provides the driving terms for static fluctuations rj s and r) v induced by 
the isoscalar and isovector electric field (the terms linear in the fluctuations 
do not vanish)Q To guarantee these fluctuations to be orthogonal on the zero 
modes z c a we have to implement the constraints 

J d 3 r z c a n 2 ab ri^ v = for all c (3.90) 

simply by adding them with Lagrange multipliers Af and A)f to the la- 
grangian. The relevant zero modes are the three infinitesimal translations 

(3.91) 

lr This was first noticed by B. Schwesinger 
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normalized appropriately ( |2.36| ) with Mo being the classical soliton mass 
( |3.3| ). Mathematically fluctuations unquestionably have to be defined in the 
space orthogonal to the one spanned by the collective coordinates (Dirac con- 
straints), physically the constraints (|3.90| ) make sure that simple translations 
of the system may not contribute to the electric polarizability. 

With these preparations we are ready to solve the e.o.m. 



h 2 ab Vb 



5Bc 



(er 



+ K n abZb 



77=0 



(3.92) 



h 2 abV v = e{er)D 



3p 



6V$ 



Si] a 



+ ^n 2 ab z c b 



(3.93) 



77=0 



for the components induced by the electric field, which reinserted into the 
lagrangian 



"77 



e -Jd 3 r(erM + ^) (d 3 




77=0 



8r] a 2 5r] a 



Tj=0 



tSS 1 t SV _i_ tVV 



SV 



(3.94) 



lead to terms proportional to e 2 which contribute to the electric polarizability. 
By looking at the Legendre transformation we made sure that this naive 
insertion into the lagrangian is here as an exception of the rule correct. 

All three terms in ( |3.94| ) are 0{N~ X ), L?r and VOf contribute to a and 
only LzY contains an angular velocity which splits neutron and proton electric 
polarizabilities. In contrast to a where we have many other contributions in 



0(N C ) (e.g. 2-loops) the leading contribution to a« is contained in ([3.94 ) 



-sv 



eD 



3p 



Sr 



er 



5V Q p 



Si] a 



Va 



e f jo , ,8B 

2J dr(£r) Wa 



T7=0 

V 
Va ■ 



(3.95) 



77=0 



Although we have calculated both isoscalar and isovector induced compo- 
nents and checked our calculation using the equivalence of the two expres- 
sions in ( |3.95| ) which follows from the e.o.m., we shall present here only the 
evaluation via the first one using the isoscalar component in some detail. 
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Figure 3.4: Radial functions u(r) and v(r) 
of the induced isoscalar electric dipole mode 
plotted over a range of 3/m. The node of v(r) 
is due to orthogonality on the translation. 



With Bq expanded linearly in the fluctations, the e.o.m. for the induced 
isoscalar component follows from ( [3.92 ) : 



— [e„ - r a (re)) + %f a (re) 



+ ^n 2 nh z c h . 



c lb ab^b 



(3.96) 



The Lagrange multipliers Af = \ee c / y/Mo, determined by multiplying ( p.96| ) 
with zz from the left are inserted into the e.o.m. 



hlbVb 



e 



+ 



f 2 

J IT T^l 



2vr 2 r 2 ' M 



+ ±Z-F'b L )r a (re) 



(3.97) 



From this equation it is noticed that the electric field induces in this case a 
pure El mode 



[u(r)r(re) + v(r)(e — f(fe))} 



(3.98) 



47T 2 / 2 

and it is evident that it is the translational zero mode which enters the 
constraints Q3.90| ) . The radial functions u{r) and v(r), subject to boundary 
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conditions u'(0) = v'(0) = , u(oo) = v (oo) = 0, are depicted in fig. |3^| v{r) 
has a node which is caused by orthogonality on the translational zero mode. 

The El mode i] s yields after averaging over the electric field directions 
(first expression in ( |3.95|) ) 



2e 2 f°° . 3 



9n z f n i3 Jo 



2 s 3 

sc b T u + (eg - 4c4 + c 6 F' )— (cw - v) 



-c^ 3 m + (cl - 2c\ + c^)F's 2 u' 



(3.99) 



Here we used that the isospin L 3 = —QD 3a Q^ is related to the angular 
momentum Q^. 

After presenting the mechanism which splits neutron and proton electric 
polarizabilities with a positive a« of 0{N^< 1 ) we will give an analytical result 
in the chiral limit — > 0. It turns out that aoj diverges in this limit such 
that the integral in ( |3.99| ) becomes dominated by the asymptotic region for 
which we are able to solve the differential equations ( |3.97|) analytically 

u = — — — [2 + m-r)e , v = — — — e " . 3.100 

With these solutions we obtain 

m^o 2e 2 [°° 3 e 2 g\ A 

m ~* o 2f o / rfr r Fu = T"^T~^ XT' 3.101 

971-^6 Jo 47r Qnf^m^ M 

where A represents the nucleon-A split. Physical values for f n , Mo, A 
and m n result in a« = 12.5 • 10~ 4 /m 3 , which is too large. However, the ratio 

m m ^ 16 A 1 

a 15 M 3 V ; 

turns out to be not unreasonable. Numerical values for the exact expression 
(|3.99|) are given in the results chapter. 



3.5.4 Magnetic polarizability 



Finally we will present a complete calculation of the static magnetic polariz- 
ability (3 to leading order O(Nq). All previous attempts plf |2"3fl |24]| miss an 
important part which may be noticed by the violation of the HBChPT re- 
sult (3 " = a/10 @@.In the course of our derivation we will recover this 
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relation although unfortunately in the end we will loose it again for reasons 
which will become obvious. Consequently the magnetic polarizability is the 
only quantity where the HBChPT result in the chiral limit is not contained 
in our tree approximation. 

There is agreement on the contribution due to the seagull term (see dis- 
cussion in section 3.5.1) which is the term quadratic in the magnetic field b 
when the isovector photon field vf — ~(r X b)i is inserted into the lagrangian 



2„2 



fie 



120 



D 3a D 3b I d\r 2 s 2 



(b T -ct-)(6b 2 S ab + 2b a b b ) (3.103) 



+ (ct-2cl)-(b 2 5 ab + 7b a b b 



vy 



(2c\ - c{) J d\ s 2 b 2 . 



Using D 3a D 3b 
bution 



^5 ab for nucleon states we obtain the direct 0{Nq) contri- 



py 



d/ir 



367T 

f_y 

3tt 



dr r s 



2„2 



(2c* - c{) J d 3 r s 



2c 



(3.104) 



to the magnetic polarizability. The last term again represents the non- 



minimal couplings which vanish for c 4 



2c\. 



In the chiral limit we obtain 



Pdir m == + —a/2 (compare ( [3.71| )) which is not the desired HBChPT result 
indicating a missing contribution. 

In analogy to ( |3.89|) we consider the term linear in the magnetic field 
(non-minimal couplings do not contribute) 



:, = -J d 3 r(rxb) l (D 3a V l a + -B l ). 



(3.105) 



Because the magnetic field couples via the spatial components of the vector 
current, the isovector part is leading O(N^) in contrast to the electric field 
case in ( |3.89| ) (but similar to the coupling of the electric field in vector meson 
lagrangians) and will therefore give a contribution to the magnetic polariz- 
ability in leading order. As in the electric case of the previous subsection the 
isoscalar part of O^Nq 1 ) splits the polarizabilities of neutron and proton. In 
the following we concentrate on the leading O(N^) term. 

Previous attempts consider this term with the spherical hedgehog inserted 



3e 
2M 



\i D 3p b 



3p u p 5 



(3.106) 
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in second order perturbation theory which leads to the standard paramag- 
netic A contribution 



Pi 



para 



e 2 (H V ) 2 
M A - M N 



(3.107) 



This procedure overestimates the effect of Q3.105T ) considerably (see table 
|3.2| ) because it neglects the soliton's response to the external magnetic field. 
Although /3 para has the correct 0(N}j) it remains finite in the chiral limit and 
cannot add to the relation (3 miT = a —a/2. Therefore it is not the desired term 
we are looking for. 



The magnetic field in ( |3.105| ) gives rise to deformations of the soliton. 
Because we are able to make the magnetic field arbitrarily small, it suffices to 
calculate the linear response via the vector current in ( |3.105| ) which provides 
the driving term for the static fluctuation 



KbVb 



~(r x b) t D 3p — 



+ ^c n a\b Z b 



(3.108) 



77=0 



D 3p [2scr a (b p - f p (fb)) - 2sf p (b a - r a (rb)) + 



+ ^ n2 a b Z b ■ 

The N£a model terms are listed explicitly and the zero mode constraints 
(compare previous subsection) are added for later convenience. The solution 
reinserted into the lagrangian leads 



-vv 



Sr r} v a 



+ ~ 2 (r X b),D 3p 



r 7=° 



(3.109) 



eft 



D 3p J Sr [2scr] L (b p - f p (fb)) - 2sf p (br7 a 



-Uw(3 ind )b 2 



to a term proportional to b . Lets forget about the constraint for a moment 
and solve ( |3.108| ) analytically in the asymptotical region. The result 



V r— >oo 

Va = —(^p> e D 3p{rpb a - r a b p ) 



(3.110) 



inserted into fl3.109| ) leads after some simple integrations to f3 — 3a/5 

which is exactly the HBChPT result (3 = (3dir + fiind m ^*° a/10 as promised. 



However the full equation (|3.108|) cannot be solved without the constraint 
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because the rotational zero mode is solution to the homogeneous equation. 
More explicitly this is seen by multiplying ( |3.108| ) with z% from the left 



K = -f V^VqA , (3-111) 

which determines the lagrangian mulipiers 7^ 0. It is worth noting that 
fl3.11ip does not only hold for the N£a model but also for the full lagrangian. 
Apart from the mathematical necessity of the constraints in eq. ( p.!08| ) they 
are also needed for physical reasons: we do not want simple rotations of 
the soliton contribute to the magnetic polarizability. When the lagrangian 
multipliers ( |3.111| ) are inserted into the e.o.m. ( |3.108[ ) 

h lbVb = — ^D 3p [2scf a (b p -f p (rb))-s(5 ap (fb) + b a f p -2f a f p (fb)) 

+ ■••], (3.112) 

it is noticed that the magnetic field induces a monopole and quadrupole 
deformation 

vl = -^ r D 3p [2f(r)b p + u(r)(b p -3r p (fb))] (3.113) 
VaT = -7T D 3 P v(r) [S ap (rb) + f p b a - 2f a f p (rb)] , 

while the magnetic dipole is projected completely from the e.o.m. ( |3.108| ) 
by the constraint. The radial functions subject to the boundary conditions 
/(0) = u(0) = v (0) = and 

jyr) = u{r) = v(r) = e w (3.114) 

(compare ( |3.11L)p ) are depicted in fig. ^75] Because of the constraint we now 
obtain [3i n d rn " T =° 9a/20 and (3 —a/20 in the chiral limit in variance with 
the HBChPT result j55j 6q| . We could not decide whether this discrepancy 



is located in HBChPT which does not know zero modes, or is caused by the 
soliton model's rigid rotator approach which rotates the profile rigidly for all 
distances even for r — » oo. 



Table |3.2| lists the complete leading 0(Nq) contributions for the param- 



eter sets A, B, and C. For comparison we also give the perturbation theory 
result P P ara (|3.107| ) which overestimates j3 in d considerably and may at best 



be considered a rough estimate. From this table it is noticed that the total 
magnetic polarizability comes out small in magnitude for all three parameter 
sets due to a strong cancellation of the direct and induced contributions. Due 
to this cancellation in leading 0(Nq), 1-loop corrections of 0(Nq) may play 
a decisive role for this quantity. 
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Figure 3.5: Radial functions /(r) of the 
monopole deformation and u(r) and v(r) of 
the quadrupole deformation induced by a con- 
stant magnetic field (full lines). For compari- 
son the asymptotical function relevant in the 
chiral limit is also plotted (dashed line). 



However loop corrections for the magnetic polarizabilities are not pre- 
sented in this paper for the following reasons: even if we were able to cal- 
culate the 1-loop contributions to fi^ir (|3.1Q4j) , where we expect a drastic 
reduction of the absolute tree value similar to the electric polarizability in 
subsection 3.5.2, we are not yet in a position to calcuate loop corrections to 
a deformed soliton necessary for [3 ind . Because both loop corrections are of 
the same O(N^), it does not make sense to consider only one of them alone, 
still it is likely that the loop corrections will shift the magnetic polarizability 
in table |3.2| to small positive values. 

Because experimental information on static polarizabilities is scarce, it 
would be appropriate to calculate directly the Compton scattering ampli- 
tudes within chiral soliton models in order to avoid the model dependent 
and uncertain extraction of the static polarizabilities from these data. In 
principle this should be possible. 
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Table 3.2: Direct and induced contributions 
to the magnetic polarizabilities are compared 
for the three parameter combinations A, B, C. 
The perturbation theory result I3 para overesti- 
mates f3i n d by a large margin. All polarizabili- 
ties are given in units 10 _4 /m 3 . 





A 


B 


C 


ftdir 


-8.6 


-9.4 


-10.7 


Pind (EO) 


2.7 


2.5 


2.0 


(3 i n d (E2) 


5.4 


5.9 


6.8 


ftind 


8.1 


8.4 


8.8 


Ppara 


11.0 


15.1 


23.4 


P 


-0.5 


-1.0 


-1.9 



3.6 Electromagnetic properties of the A iso- 
bar 

To complete our discussion of electromagnetic properties we briefly adress 
some of these quantities connected with the A isobar. In general the static 
properties of the A are related to those of the nucleon by simple geometrical 
coefficients which follow from the collective operators involved. For adiabatic 
quantities this operator structure remains unaffected by 1-loop corrections 
and consequently the model independent relations survive for these quan- 
tities. The experimental verifications of these relations may be considered 
as a crucial test of the validity of the l/Nc expansion. As an example we 
discuss the magnetic moments in the following paragraph. There we also 
mention the quadrupole moment (which vanishes for the nucleon) as well as 
the corresponding transition moments. 

Another interesting quantity subject to recent measurements is discussed 
in the second subsection, namely the electromagnetic ratio of the A photo- 
decay amplitudes. All soliton model attempts to calculate this quantity so far 
make use of Siegert's theorem at the finite photon momentum where it is not 
stictly valid. We show that the correct treatment reduces the electromagnetic 
ratio considerably. 
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3.6.1 Magnetic and quadrupole moments 



The operators for the calculation of magnetic and quadrupole moments are 
given by 

fi = -3fi v D 33 ~2fi S R 3 + 0(Nc 1 ) 

Q = -\<r 2 > e (D 33 R 3 -^D 3p R p ) , (3.115) 

with the isoscalar and isovector magnetic moments fi s and fi v and the 
squared radius < r 2 >e of the moment of inertia's density defined in section 
3.4. The adiabatic quantity fi V may contain the 1-loop corrections described 
in the section just mentioned. Evaluating the matrix elements of fj, ( [3.115] ) 
with respect to A states or nucleon and A states 

^ ++ = <A+\S 3 = ^\A+\S 3 = l>= 9 -v v + 3v s + 0(Nc 1 ) 

^ NA = <A + ,S 3 = ^\p,S 3 = ^>=V2fi V + 0(Nc 1 ), (3.116) 

relates the magnetic (transition) moments by simple geometrical coefficients 
to the isoscalar and isovector magnetic moments. Upon using the experimen- 
tal values fi v = 2.35 and fi s = .44 we obtain fi A++ = 5.55 and fi NA = 3.32. 
If accurate experimental information about these magnetic moments were 
available we had a check on the 1/Nc expansion i.e. we could estimate the 
next-to-next to leading order contribution generated by the vector current. 
In our model such 0(N^ 1 ) contributions are produced in tree approximation 
by small soliton deformations proportional to the angular velocity squared 
(next subsection), by nonadiabatic 1-loop processes which come with one an- 
gular velocity, and by 2-loop diagrams (compare a similar discussion for the 
axial current in appendix A). Unfortunately there are no precise experimen- 
tal data: for /i A++ the particle data group gives a range of 3.7 — 7.5 nuclear 



magnetons ||78|| , for the transition magnetic moment there exist photopro- 
duction data, but these may not directly be compared to the matrix element 
( pT6|) (see table 3.4). 

Many other A observables like rms radii and polarizabilities are related 
to those of the nucleon by simple geometrical coefficients just as in the case 
of the magnetic moments. We do not list them here because there are no 
experimental data available, instead we concentrate on the electromagnetic 
ratio of the A photo-decay amplitudes. For this purpose we add here still the 
quadrupole moments 

Q A++ = <A+\S 3 = 3 -\Q\A+\S 3 = 3 ->=~<r*> e 
Q NA = < A\S 3 = ±\Q\p,S 3 = ±>=-^<r 2 > e . (3.117) 



Our favorite models A and B give Q A++ — — 0.06/m 2 , the transition mo- 
ments are listed in table 3.3. Note in this context that the quadrupole mo- 
ments are nonadiabatic quantities of 0(Nq) which we calculate in tree ap- 
proximation only. 



3.6.2 E2/M1 ratio of the A photo-decay amplitudes 

The desire for precise values of the A photo-decay amplitudes has spawned 
numerous theoretical and experimental investigations in recent years. Much 
of this effort has been directed toward the E2/M1 ratio of these amplitudes. 
In soliton models from the earliest attempts on |f79| ||80|| ||8l| down to very 



recent publications |B2|p3[ the E2 transition was always calculated via the 



time component of the vector current Vq by employing Siegert's theorem |84 



at the finite photon momentum, where it is not strictly valid. We calculate 
the E2 transition directly from the spatial components of the vector current 
V?, show that indeed for q — > the transition matrix element may be ex- 
pressed by the corresponding matrix element of Vq (Siegert's theorem), but 
we find nonnegligible corrections at the photon point leading to a noticeable 
reduction of the electromagnetic ratio. 

For this purpose we have to solve the e.o.m. 



- d% a = dtV? = V a = ^s 2 b T D ap (f x R) p (fR) (3.118) 

to O^N^ 1 ) because the E2 transition is supressed by 1/Nq as compared to 
the Ml transition. This leads to small rotationally induced monopole and 
quadrupole deformations 



Vl 



2f(r)R 2 + u(r)(R 2 -3(fRY) (3.119) 



1 



= -^j^v(r)(R-r(fR))(rR) 

described in appendix A.l. The system of differential equations for the ra- 
dial functions /(r), u(r) and v{r) is identical to the one discussed for the 
magnetic polarizabilities (subsection 3.5.4), and consequently also the radial 
functions shown in fig. 3.5. With these components inserted, the vector cur- 
rent obeys ( |3.118|) . For this reason we will not encounter ambiguous results 



calculating the E2 transition via the time or the spatial components of the 



electromagnetic current, discussed for constituent quark models in |85[ . 

In order to minimize errors, the so-called "equal- velocity" frame where 
the incoming nucleon and the outgoing A have opposite velocities and which 
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reduces to the Breit frame as the nucleon-A split vanishes, is chosen as a 
convenient reference frame [86]. In this frame the photon momentum is q& = 
(M A - M N )/2^M A M N = 296MeV. We define the electromagnetic ratio by 
the ratio of the corresponding helicity amplitudes 



E2 lAi(E2) i<A+,S a = ±\M&\p,S a 



Ml 3Ai(Ml) 3<A+,S 3 = ±\MM\\p,S 3 



(3.120) 



with the multipole operators |87 

M A M1 (g 2 ) = iVQnX I d\V? h{qr) 



M\ (q ) 



VIO 



7T 



q 



J d\V? [V x (j 2 (gr)r 22A (f))] 4 . (3.12i; 



First we are going to evaluate the Ml transition operator by inserting the 
vector current in tree + 1-loop (which part may be expressed by the isovec- 
tor magnetic formfactor discussed in section 3.4) as well as the rotationally 
induced soliton deformations of 0(Nq 1 ) 



3A 
2 



M N 



Srji- (c(5 pX - f p f x )r]L - r p r] TX H ) 

r 



3A 


q G 


T 


[m n 




1 




45@2 




1 


' 150 2 



'3A 



(3.122) 



/ d\-j 1 ~(2cu + 3v + ---)L 3 R > 
J r 



With this multipole operator, the Ml transition matrix element in tree + 
1-loop with soliton deformations included 



2 

Mv 



> 



q 



(3.123) 



G%W) = V2{G V M (q 2 ) + 20iie2 



is determined. 



2V2M N 

J d 3 r 3l S -[2c(10f -u)-3u + •-•]} 



For the E2 transition there is no choice: the soliton deformations ( |3.119| ) 
must be inserted, otherwise the operator vanishes. Alternatively, and of 
course with the same result, we may employ partial integration and the vector 
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current conservation fl3.118| ) 




2\ 



(3j 2 - qrj 3 )V* (3.124) 



+f 7T q 2 j2D 3p (rF'cr x 77 — rsf x 77' 



2X 




q 2 



f 2 (3j 2 - qrj 3 )s 2 b T - ^j 2 (rF'cv - rsv' 

D 3p (R p - r p (fR))Y 2X . 

In the last step the soliton deformations ( p.H9|) (only N£a contributions are 
written out in explicit form) are inserted, and it is then noticed that the 
entire operator may be written as a total time derivative which in the end 
may be replaced by iq (compare the derivation of Siegert's theorem). The E2 
transition matrix element is now 

< A+,S 3 = 1 -\M^\P,S 3 = -\ >= -|^G^(g 2 ) (3.125) 

y/2 - r ~ 2 



f 2 (3]2 - qrj 3 )s 2 b T —(rF'cv - rsv' + 



q 32, 

9q 2 Q 

It should be remembered here that is a nonadiabatic quantity which 

we calculate in tree approximation only. There may possibly be nonadiabatic 
1-loop contributions of the same order which to calculate we are not yet in 
a position. 

For q — ► we obtain — > Q NA ( p. 11 7[ ) which is indeed Siegert's theo- 
rem. At the photon point the corrections (second and third term in (|3 . 1 25|) ) 
are sizeable, numerical values for the parameter combinations A, B, and C 
are compared in table 3.3. The reason for \G^ 2 ~ Q I being larger than 
expected is connected with the fact that the induced rotational components 
(fig. 3.5) fall off quite slowly which corresponds to a large effective radius 
of these corrections. With ( [3.123j , |3.125|) the electromagnetic ratio Q3.120 ) is 



determined as 

E2 q A MG^(ql) q A M G^(A> 



Ml 2 G^(ql) 2 G«fi(A*) 



(3.126) 



In the last step a relativistic correction similar to (3.32) leads to the replace- 
ment q A — ► A = Ma — = 293MeV in the transition formfactors which 
numerically is of negligible influence. The results are listed in table 3.3 with 
an absolute value considerably smaller than in previous soliton calculations 
mainly for two reasons: 
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Gmi is larger and close to the experimental value at least for our favored 
models A and B (accidentally, maodel A hits the experimental value of 
the corresponding helicity amplitude |78| exactly) because we included 



loop corrections of 0{N^) and soliton deformations of 0(N C 1 ] 

• The consideration of G^(A 2 ) instead of Q NA (A 2 ) (deviation from 
Siegert's theorem due to finite A) reduces the absolute value of the 
electromagnetic ratio further. 

Our result of E2/M1 ~ —1.7% has to be compared with the experimental 
value reported by the particle data group — 1.5 ± 0.4% and a recent 
MAMI experiment with the preliminary result — 2.4 ± 0.2% 



Table 3.3: Magnetic and quadru- 
pole transition moments in units 
of nuclear magnetons and fm 2 
resp. are compared for the three 
parameter combinations A, B, 
C. The corresponding transition 
formfactors are evaluated at the 
photon point. 





A 


B 


C 


/i ivA 


4.1 


4.4 


4.9 




3.2 


3.4 


3.8 


Q NA 


-.034 


-.035 


-.038 


Q NA (A 2 ) 


-.021 


-.022 


-.023 


G£ 2 A (A 2 ) 


-.015 


-.015 


-.016 


mm 


-1.7 


-1.6 


-1.5 



Because the extraction of the electromagnetic ratio from experimental 
data is far from being trivial [39] [3C], a complete calculation of the 'yN — > ttN 
reaction, which includes the background with the A resonance lying in the 
continuum would be desirable in order to obtain complex helicity amplitudes 
for the NA transition. Although such a calculation should in principle be 
feasible within soliton models, the effort would be tremendeous. 
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Chapter 4 
Results 



In this chapter, we are finally going to present and discuss results obtained by 
applying the method outlined above to the quantities introduced in the pre- 
ceding paragraphs. As a shorthand, we shall in the following denote quantities 
whose loop corrections can be treated in adiabatic approximation (and are 
therefore calculated in this paper) as adiabatic quantities, whereas the others 
(for which we are unable to calculate the loop correction) will be summarized 
as nonadiabatic ones; only their tree level values will be given. For all quan- 
tities considered, we evaluate the Lagrangian ( |2.47| ) for the three parameter 
sets A (e = 4.25, g u = 0), B (e = 4.5, g u = 1.0) and C (e = 5.845, g u = 2.2) 
introduced in section 2.2 . From the remarks in that section, it should be 
clear that we cannot expect model C to result in sensible numbers since the 
local approximation with vector meson parameters overestimates the effect 
from the u) meson and neglects important contributions of higher chiral or- 
ders. However, sets A and B, designed to simulate these missing pieces by 
introducing 'effective' parameter values should give an accurate picture of the 
abilities (and defects) of the soliton approach with respect to the quantities 
under investigation. Since amongst these three, set B comes closest to the 
original vector meson model, we consider this to be the most realistic one. 

Before entering the discussion of magnitude and sign of the 1 loop correc- 
tions it is worthwhile to recall the salient features of the tree approximation 
within the three parameter sets A, B and C, at least for the adiabatic quan- 
tities: 



The soliton mass has been a constant source of embarrassment for 
the past decade, since it came out too large by a factor of more than 
3/2 almost regardless of parameter combinations used or meson species 
included in the model. This is also obvious from tables 14. 1L 14.21 and [O 



where the tree level mass is shown to be in the 1500 — 1630MeV range 
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for the three models. 



With values between 54 and 67MeV against experimental 45MeV, the 
<t term is also overestimated in tree for all parameter combinations 
under consideration here. This is not a common feature of all Skyrme 
type models and depends on our usage of the full ChO 4 lagrangian 
including symmetry breakers. The effect is much less pronounced than 
in case of the mass, especially if the comparably large error margin of 
the experimental datum is taken into account. Nevertheless, a sizeable 
reduction due to 1 loop corrections would be necessary. Partly due to 
the overestimated a term, the scalar radius comes out consistently much 
too low in tree, its value of around lfm as opposed to experimental 
1.6 fm being rather insensitive to the parameters used. 

The axial vector constant is a weak point of Skyrme type models, 
since the tree contribution doggedly refuses to grow appreciably beyond 
1.0 (our values being .91 — 1.03) as long as reasonable parameters are 
considered, where reasonable means that the other quantities should 
not be completely off the mark. Likewise, the quantity g\ < r 2 >a is 
underpredicted in all cases, numerical values ranging from .36/m 2 to 
.42/m 2 whereas experiment points to .53/m 2 . 

Concerning the electromagnetic formfactors, it has long been rec- 
ognized that magnetic moments generally come out too small, most 
notably the (nonadiabatic) isoscalar one, which is wrong by a factor of 
2, but also the isovector moment for which our parameters give 1.6 — 1.8 
nuclear magnetons in contrast to the experimental finding of 2.35. Sim- 
ilarly, for some time, conventional wisdom had it that e.m. radii were 
seriously underpredicted in purely pionic models. However, inclusion 
of nonminimal couplings with coupling constants chosen in accordance 
with the logic presented in section 2.2, as well as taking into account rel- 
ativistic corrections serve to remedy this problem already at tree level, 
and for our parameter sets, actually the converse is true: isoscalar elec- 
tric as well as isovector magnetic radius are slightly overestimated with 
values of .61 — .67 fm 2 for the former and .77 — .72/m 2 for the lat- 
ter case to be compared to experimental data of .59 fm 2 and .73 fm 2 , 
respectively. 

The electric polarizability is consistently too large in tree without 
regard of the model considered; the discussion in the corresponding 
section of the previous chapter has revealed that this would also hold 
in models containing vector mesons explicitly, which were for some time 
believed to cure that deficiency. Actual numbers are in the 18 — 22 • 
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10~ 4 /m 3 bracket for our models; however, the amount by which the 
polarizability is overestimated is difficult to quantify in view of the 
large error bar, the experimental value being 9.5 ± 5 • 10~ 4 /m 3 . 



Turning now to the loop corrections, inspection of tables |4.1| , ^4.2| and [13 
provides a rough classification in terms of the sensitivity to the stabilisation 
mechanism: quantities with a loop contribution depending sensitively on the 
choice of parameters comprise a term, g^, axial radius end electric polariz- 
ability. All others are rather insensitive to the details of the lagrangian. We 
shall postpone discussion of the first kind of properties and deal firstly with 
the latter species. Here, the most pronounced results are the following: 



• As had been shown previously ( |[36| ||35|| ) the soliton mass is consider- 
ably reduced towards its experimental value once loop corrections are 
taken into account. Every parameter set considered here results in a 
tree + 1 loop mass of around 950 MeV. 

• No less importantly, the scalar radius is considerably enhanced to- 
wards the experimental datum, with tree + 1 loop results in the range 
of 1.3 - 1.5/m 2 . 

• The same is true for the isovector magnetic moment, which shows 
an enhancement to 2.24 — 3.13. The latter number (3.13 for set C) ac- 
tually reverses the tree level underprediction into a considerable over- 
estimation, but both the other models are satisfactory. 

In all these cases, the quantum correction has the right sign and a magnitude 
of slightly more than 1/3 relative to the tree level value. This serves as an 
impressive justification of the conjecture that loop corrections might play a 
crucial role in the understanding of baryon properties in soliton models. 

Less marked changes occur in the e.m. radii, where the correction is around 
1/6 of tree, yet has the right sign in all cases except C, where the too small 
tree magnetic radius is further reduced. Since the tree level values were al- 
ready rather good, it is not really justified to speak of an improvement here; 
however, the fact that loop contributions are tiny where one would expect 
them to be, and generally have, once again, the right sign supports the hy- 
pothesis that quantum corections, properly taken into account, could even- 
tually cure the defects of tree level evaluated Skyrme type models without 
worsening things with respect to quantities where tree already gives accept- 
able results. 

So far, there has been little to choose between sets A, B and C with the 
exception of magnetic properties in case C. This radically changes if one 
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considers the first kind of quantities mentioned above, namely a term, g^, 
axial radius and electric polarizability. 

Dealing first with the <x term, this happens to be corrected in the right 
direction and with the right magnitude for the first two parameter sets A 
and B; actually, set A turns the too large tree value into a too small tree + 
1 loop one (32MeV), whereas set B is just about right (44MeV). However, 
the startling observation is that for a strong enough sixth order term (set C) 
the 1-loop contribution changes sign and increases the already too large tree 
value to about 80MeV. A similar phenomenon occurs in case of the g A if one 
takes into account set D, and also for the axial radius. We do not want to 
repeat the discussion of section 3.3 here, but once again we emphasise that 
the changing sign of the loop correction serves as a filter to decide which 
model to discard and which to keep; sets C and D are clearly ruled out by 
the numbers. (Note that in contrast to the situation for the slope of the 
axial formfactor g& < r 2 >a, the 1 loop correction to the radius itself is 
not worrisome, and that < r 2 >a comes out too large, yet within the error 
margin of the experimental datum, if one includes the estimated 1/Nq piece). 
The fairly good stability of the tree + 1 loop value of gA over a wide range 
of scale as shown in fig. |3.3| probably rules out missing higher chiral orders 
to result in a scenario different from the one presented here (e.g., positive 
loop correction for otherwise reasonable parameter sets). Eventually, one of 
the success storys of our endeavour is the loop contribution to the electric 
polarizability, which turns out to bring both, models A and B, to the 
experimental datum including error bars (tree + 1 loop values being 9.8 and 
14.1 respectively). Like its precursors discussed above, the correction changes 
sign as one proceeds from set B to set C, inflating the tree value of 22 to 
about 27, therefore adding further to the list of its deficiencies. We believe 
this sensitivity to the strength of the sixth order term in the lagrangian to 
reflect the problems of the local approximation on the uj meson discussed 
above: A strong sixth order term gives rise to unphysical features in the 
scattering phaseshifts, and consequently fails to result in sensible predictions 
for certain quantities. Since we cannot, for the moment, properly account for 
the uj meson, we have to be content with the philosophy adopted in section 
2.2, namely to use effective values for e and g u which then results in good 
agreement with experiment. 

Summarising, we have shown loop corrections to bring all adiabatic quan- 
tities with the exception of the axial ones, close to their experimental values 
for parameter sets A and B. In case of g A and < r 2 >a, we believe (section 
3.3) that their status is exceptional and that the 1/Nq expansion must run 
into difficulties in their case. We have provided an estimate of the next- to- 
next-to-leading order correction to them, which turned out to be large and 
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Table 4.1: Tree and 1-loop contribution to various quan- 
tities for parameter set A (e = 4.25, = 0). 





tree 


1-loop 


E 


exp. 


M[MeV] 


1629 


-683 


946.0 


939 


a [MeV] 


54 


-22 


32 


45 ± 7 


<r 2 > s [fm 2 ] 


1.0 


+.3 


1.3 


1.6 ± .3 


9a 


.91 


-.25 


.66 


1.26 


<r 2 > A [fm 2 ] 


.45 


-.04 


.41 


42 +- 18 

-.08 


<r 2 >E if™ 2 } 


.62 


-.11 


.51 


.59 




1.62 


+.62 


2.24 


2.35 


<r 2 > v M [fm 2 ] 


.77 


-.13 


.64 


.73 


a [10~ 4 /m 3 ] 


17.8 


-8.0 


9.8 


9.5 ± 5 



Table 4.2: Tree and 1-loop contribution to various quan- 
tities for parameter set B (e = 4.5, g w = 1.0). 



M[MeV] 



a [ MeV] 
s [fm 2 



<r 2 > 



9A 



<r 2 > A [fm/ 



<r 2 > s 



/i 



v 



<r > M 



[fm 2 



a [10 A fm c 



tree 



1599 



58 



1.1 



.96 



.44 



.64 



1.69 



.76 



19.4 



1-loop 



-646 



-14 



+.3 



-.15 



-.01 



-.09 



+.1 



-.17 



-5.3 



E 



953 



44 



1.4 



.43 



.55 



2.57 



.59 



14.1 



exp. 



939 



45 ± 7 



1.6 ± .3 



1.26 



.42 



.08 



.59 



2.35 



.73 



9.5 ± 5 



Table 4.3: Tree and 1-loop contribution to various quan- 
tities for parameter set C (e = 5.845, g w = 2.2). 



M[MeV] 



a [MeV] 
* [fm 2 ] 



<r 2 > 



9A 



<r >a [fm 2 



< r >e [f m ' 



<r 2 > 



M [f 



m 



a [10" 4 /w 



tree 



1490 



67 



1.1 



1.03 



.40 



.68 



1.75 



.72 
21.8 



1-loop 



-539 



+ 13 



+.4 



-.11 



-.16 



-.06 



+1.38 



-.19 
+5.4 



E 



951 



1.5 



.92 



.56 



.62 



3.13 



.53 
27.2 



exp. 



939 



45 + 7 



1.6 ± .3 



.42 



1.26 

TTTST 

-.08 



.59 



2.35 



.73 
9.5 ± 5 
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remedies the defects that show up including only 0(1). Although there is in 
principle the possibility of sizeable O^N^ 1 ) effects for the other quantities, 
too, already at this point, we can safely state that whatever parameter set (if 
any) might come out of a consideration including these must lie in the vicin- 
ity of A or B: C-like sets are clearly ruled out from the electric polarizability 
alone since they would require a 1/Nq correction of around —10 as opposed 
to +5 at 0(1). Table [O gives tree results for nonadiabatic quantities. Here, 



Table 4.4: Comparative listing of tree + 1 loop 
values for all adiabatic quantities considered for the 
three parameter combinations A, B, C. Axial quan- 
tities include estimated 1/N C piece. 





A 


B 


C 


exp. 


M[MeV] 


946 


953 


951 


939 


a [MeV] 


32 


44 


80 


45 ± 7 


<r 2 > s [fm 2 ] 


1.3 


1.4 


1.5 


1.6 ± .3 


9a 


1.20 


1.29 


1.36 


1.26 


<r 2 > A [fm 2 




.55 


.54 


.59 


42 +J8 


<r 2 >E [fm 2 




.51 


.55 


.62 


.59 




2.24 


2.57 


3.13 


2.35 


<r 2 > v M [fm 2 ] 


.64 


.59 


.53 


.73 


a [lO'^frn 3 ' 




9.8 


14.1 


27.2 


9.5 ± 5 



Table 4.5: Comparative listing of tree values 
for all nonadiabatic quantities considered for 
the three parameter combinations A, B, C . 





A 


B 


C 


exp. 


A [MeV] 


290 


278 


265 


293 


M np [MeV] 


2.0 


1.7 


1.4 


2.05 ± .30 


<r 2 >E [fm 2 ] 


.86 


.90 


.98 


.85 ± .03 


n s 


.18 


.19 


.21 


.44 


<r 2 > b M [fm 2 




.74 


.75 


.77 


.59 ± .02 


A« [10~ 4 /m 3 ' 




1.7 


1.8 


1.9 


5.0 ± 5.0 



one observes that all quantities except two are rather well reproduced already 
at this level. The two exceptions are the isoscalar magnetic moment and the 
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neutron-proton-split of the electric polarizability. The latter is less severe in 
view of the huge error of the experimental value, but fj,g is seriously defective 
and it seems unlikely that a next order quantum correction would be able 
to remedy this discrepancy. In principle, however, there is the possibility of 
a quantum correction of the same order as the tree level in 1/Nq since the 
time derivative contained in the spatial component of the winding number 
current expanded to quadratic order in the fluctuations may act on a fluctu- 
ation instead of a collective coordinate and twin fluctuations count as 1/Nc 
like the angular velocity. We are not yet in a position to decide whether this 
is the root cause for the too small isoscalar magnetic moment. 

Finally, we shall briefly recapitulate the course of this report: In the be- 
ginning we posed the question whether it is possible to renormalize a chiral 
lagrangian evaluated in a nontrivial topological sector in analogy to the meth- 
ods known as chiral perturbatuion theory for the vacuum sector. The main 
problem is the restriction to a finite number of gradients contained in the 
lagrangian, which is not a priori justifiable in the presence of a soliton back- 
ground. Upon considering a lagrangian incorporating resonance degrees of 
freedom we tried to better understand this issue and concluded that using 
the 4th order chiral lagrangian augmented by a standard 6th order extension 
including nonminimal couplings would be sufficient provided a shift away 
from the experimental values of LECs would be allowed for. Later on, we 
could add substance to this conclusion by considering the scale dependence 
of the tree + 1-loop soliton mass and the axial coupling constant Qa whose 
values were shown to be remarkably stable over a wide range of the scale. 

Starting from this lagrangian, we devised a method to calculate the self 
energy of the soliton in the presence of external fields in adiabatic approxima- 
tion. In the course of the actual calculation, we had to clear several obstacles 
only indirectly related to our subject: We clarified the treatment of zero 
modes in the presence of external fields, showed how to properly compute 
the tree value of the electric polarizability and calculated the neutron-proton- 
split of this quantity. In the end questions remain as to the convergence of 
the 1/Nq expansion for axial quantities and to the size of the next-to-next- 
to-leading order contributions. However, we once again want to stress that 
in 7 out of 9 cases, the 1 loop correction has the right sign and magnitude 
to compensate for the discrepancy between tree level and experiment, which 
is truly remarkable since there are only two parameters (e and gj) involved 
(which moreover have not been tuned to achieve some specific result). 

The fact that sets A and B designed to simulate the higher chiral order 
pieces contained in a vector meson model indeed give much better results 
than set C strongly supports models which include vector mesons explicitly. 
We decided to use a pseudoscalar lagrangian only because of the arguments 



79 



presented in section 2.3 . Nevertheless, we expect that, concerning 1-loop 
corrections, the results for a lagrangian with explicit g and uj -mesons would 
come close to those of our model B. In particular the problems with the 
axial quantities discussed in section 3.3 will probably not be cured by the 
introduction of vector mesons. 
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Appendix A 

1/Nq contribution to gj± 



Although the relation (|3.45|) is model independent, the positivity of the re- 
mainder R is not a direct implication of the charge algebra. Rather, this 
follows either from the experimental value of g& (which would be a circu- 
lar argument in our context) or from the fact that R can be expressed as a 
difference of pion nucleon cross sections which is experimentally known to 
be positive. In soliton models, it is hard to see whether R is indeed positive 
independently of the parameter set chosen. Provided this were the case and 
the 1 / Nq expansion would converge for R, then would necessarily have to 
be larger than 1 as soon as the charge algebra relation ( p.44|) were fulfilled. 
We could then use a different strategy to obtain an estimate for the l/Nc 
piece of gA, which in the end is equivalent to the naive addition of the CA " 1" 
to g\ used in section 3.3.3 .Such an estimate is the subject of this appendix. 

Soliton models obey (|3.44j ) provided nonadiabatic fluctuations which come 
with a collective angular velocity are taken into account |)3|. Relative to 
the adiabatic fluctuations employed for the 1 loop contributions, these are 
suppressed by a factor of l/Nc- As a consequence current algebra is reflected 
in gA not before 0{N^ 1 ) which, of course, is also immediately clear from 

(EH). 

Quite generally, the (spatially integrated) axial current up to 0{N^, 1 ) may 
be written as 

A" = -\gA l+0) D ai + \a^{D ai , R 2 } + b^L a R t + c^D m (A.l) 

where the Nc orders are indicated by the superscripts. 

In this object, the tree contribution proportional to the angular velocity 
squared of O^N^ 1 ) due to the soliton's rotation can be shown to be negligi- 
ble (subsection A.l). Adiabatic 2-loop diagrams cannot produce the operator 
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structure of ( |A.l ). There remain the nonadiabatic 1 loop processes: We be- 



lieve that these should contribute sizeably in 0(N C ) to eq. ( |A.1| ) because 



exactly they are responsible for the restoration of the charge algebra com- 
mutation relation ([3.44 ). However, they are too complicated to investigate 
in detail. 

Therefore, in the second subsection of this appendix, in order to circum- 
vent this obstacle, we are going to estimate the l/N c contribution to qa 
produced by nonadiabatic 1 loop processes by imposing the S77(4) current 
algebra commutation relations 



["4-3 > "4.3] — ^ ^abcL c 



(A.2) 



A.l Contribution due to rotationally induced 
soliton deformations 



Here we show that the effect of soliton deformations on the axial current is 
too weak to explain the relatively large 1/Nc contributions needed in order to 
fulfill the CA commutation relation (A.2). This does of course not imply that 
soliton deformations are small by themselves, in the vector current they are 
responsible for sizeable contributions to the magnetic moments as we have 
noticed in subsection 3.6.2. Small deformations due to the soliton's rotation 
are driven by the term 



Ln = U J d 6 r [sc(r x Q, K frfl - s (rfi*)(fi*T^ 
in the lagrangian. The corresponding e.o.m. 

may again be solved analytically in the asymptotical region 

~r(n R ) 2 - n R (fn R 



hl^ = f n \scr a (r X n R ) 2 - s(Q a - r a (rn R ))(r£l R ) + 



16nf n 



(A.3) 



(A.4) 



(A.5) 



which is in accordance with the result in ref. [91|. The full equations induce 
a monopole and a quadrupole deformation 



Vl 



6/. 



2f(r)(n R ) 2 + u(r)m R f-3(fn R f) 



R\2 



\R\2\ 



(A.6) 



l —v(r)(Sl R - r(fQ R )(rn R )) 
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and it turns out that the system of differential equations for the radial func- 
tions f(r), u(r) and v(r) is identical to the ones discussed for the magnetic 
polarizabilities (subsection 3.5.4) and the electromagnetic ratio of the photo- 
decay amplitudes (subsection 3.7.2). The radial functions are depicted in 
fig. |3.5| . These rotationally induced components may give rise to 1/Nq contri- 
butions to the spatial axial current 



A\ = J Sr x 4 (ig - d^AD 



f 2 



20 2 

2 

m 
96 

2 

m 



156 



J d\ r scb T (^{D ai ,R 2 } - L a Ri) 
| J d\rcf l -{D au R 2 } 

L J d \ r {\ cu + ±v)(~{D ai , R 2 } - 3L aJ R 4 ) . (A.7) 



Here we used partial integration and we listed only the contribution of the 
familiar symmetry breaking mass term. Equivalently, we could have inserted 
( |A.6| ) directly into A® with the same result. 

From ( |A.7|) , we recover the operator structure ( |A.1| ) and we may read off 
the tree level pieces of the coefficients a* -1 ) and Evaluating ( |A.7|) in 

the chiral limit with the appropriate asymptotic expressions Q3.114] ) 



9 A l rn „ 21 



@2 m 2 2 



{D at ,R 2 } 



9 A A f r , „ 2 



-{-Dai, R } — 3L a i?j 



26 2 m2 v 2 

= • 0{m- 2 ) (A.8) 

the divergent terms cancel exactly which is pleasing because we expect the 
coefficients c^ -1 ) and to be finite in the chiral limit in order to fulfill 
the CA commutation relation( |A.2|) . According to this cancellation, the con- 
tributions c^ -1 - 1 and £>( -1 ) remain small also for finite pion mass. For example, 
model A provides the tree contributions a^" 1 ' = — .04 and b^" 1 ^ = —.09 which 
are negligible compared to the values a^ 1 ^ = —b^ 1 ^ = .36 necessary to fulfill 
the CA relation ( |A.2| ) (compare also the following subsection). Therefore we 
do not expect soliton deformations to appreciably contribute to this relation. 



A. 2 Contribution due to nonadiabatic loops 



In this section we estimate the 1/Nc contribution to the axial current due to 
nonadiabatic loops by imposing the CA relation (A. 2). Before doing so, we 
notice that in a pure Skyrme model (N£a model with and Skyrme term 
with parameter e) the coefficient in the axial current (A.l) scale according 
to 

9\~- 2 , 9a~^- 2 , a™ ~ „ c<-i) „ _L (A.9) 
e 2 fie 2 f4 

— 1/2 

(e = 2g is 0(N C )). Thus by choosing a large Skyrme parameter we can 
make the tree + 1 loop contribution small. The property g A > 1 must then 
almost entirely be achieved due to the 1/Nc contribution. 

The CA commutation relation (A. 2) requires 

= -a(-) , (a^ - ^(g^ - V + 1 = (A.10) 

and we may simplify the expression for the axial current (A.l) 

(A.11) 



l -{D ai R 2 }-L a R l -KD a 



where we replaced C* -1 -* = — na^ 1 ^ according to the scaling behaviour (A.9). 
The constant k is parameter independent and will be determined imme- 
diately. With (A. 11) we may evaluate the weak axial nucleon-nucleon and 
nucleon-A transition amplitudes 

(1+0) , ^ (•_!) 

9a = 9 a + g Ka 1 

= g^-^ + l^. (A.12) 

With these quantities the condition (A. 10) may be rewritten as 

2 

9 A ~ ^9 a = 1 

3 -a^=g A -^g* A , (A.13) 



which is Oehme's relation [92] and a direct consequence of the CA relation 
(A. 2). From this relation it follows that qa > 1. With the coefficient a^ 1 ^ 
inserted into (A.12) we obtain 

(i+o) 4 V2 
9a = 9 A + g K {9A - -^9 a) 

4 



9a 1+0) + g<9A ~ V5i " 1) (A.14) 
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The last step involves the identification y/2/3g A = +y g\ — 1 which is rea- 
sonable because of the large Nc requirement y/2/3g A = g A . Since relation 
A.14j must hold independently of the value of g A + °^ ( c -f- the discussion in 
the beginning of this section), it follows that k > 9/4. The choice k = 9/4 
is conservative in the sense that it yields, for any given g A + °\ the smallest 
possible correction g^ 1 (fig. |A.1| ). Therefore, upon using this value for k we 
obtain 

a = i+(9.r») 2 
^ = 3 { r > , (A.i5) 

which was the estimate used in section 3.3.3 . 



Figure A.l: Plot of g A for different values of 



From the experimental value g A = 1.26 we obtain a large 1/Nq contribu- 
tion g^ 1 ^ = .40 and a small value g A °^ = .86 for tree + 1 loop. If we accept 
the assumption underlying this paragraph, namely 

• S77(4) symmetry of the spatial axial current 

then this seems to be the reason why soliton models in tree (+ 1 loop) 
approxmation always underestimate g A by a large margin. 
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